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ABSTRACT 


The  techniques  for  the  investigation  of  char  densification  phenomena 
which  had  been  developed  under  a  previous  contract  (see  BSD  TR  66-385) 
were  extended  to  include  transpiration  experiments  with  gases  and  gas 
mixtures  typical  of  the  decomposition  products  of  phenolic  resins. 

As  before,  the  chars  used  in  the  experiments  were  produced  by  charring 
high  performance  carbon  cloth  phenolic  composites  in  a  plasma  jet. 

Data  on  the  structural  characteristics  and  permeabilities  of  the  chars, 
and  the  change  occurring  during  transpiration,  are  presented.  Compari¬ 
sons  of  the  experimental  data  with  kinetic  rate  data  from  the  literature 
are  made.  A  simplified  program  for  calculating  the  change  in  permeability 
of  and  pressure  drop  across  a  growing  char  is  presented.  The  effort 
required  to  extend  this  program  to  more  complex  situations  is  discussed. 
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INTRODUCTION 


Tb*  basic  objective  of  the  study  described  in  this  report  ves  to  provide 
reliable  experimental  data  which  is  necessary  for  devising  appropriate 
analytical  expressions  for  predicting  changes  in  permeability  and  pressure 
drop  occurring  in  chars  during  reentry.  More  accurate  expressions  of  this 
type  are  needed  for  incorporation  Into  analyses  describing  transient  heat 
shield  behavior. 

In  order  to  satisfy  the  objective  of  the  study,  experimental  data  on  the 
effect  of  the  reactions  which  are  typical  of  those  occurring  in  the  char 
(at  least  for  a  phenolic  composite)  were  obtained.  These  consisted  of 
permeability  and  pressure  drop  characteristics  of  chars  as  a  function  of 
time  (time  related  to  the  reentry  time  scale).  This  data  was  correlated 
with  the  changes  in  the  internal  physical  structure  of  the  chars  and  will 
permit  better  extrapolation  of  the  results  of  the  present  project  to 
materials  of  a  somewhat  different  nature. 

Some  initial  results  have  been  reported  previously  (Reference  1).  In  that 
report  data  were  presented  on  the  changes  in  permeability  and  structural 
characteristics  which  occur  in  typical  ablative  chars  when  carbon  deposi¬ 
tion  is  caused  by  the  deconposition  of  methane  in  the  char  pore3  at  typical 
re-entry  temperatures.  For  purposes  of  con^arison,  data  were  also  presented 
for  the  densification  process  as  it  takes  place  in  commercial  graphites 
and  carbons. 

The  details  of  the  apparatus  and  techniques  developed  for  that  part  of  the 
investigation  will  be  found  in  the  aforementioned  report.  Also,  as  an 
appendix  in  that  report,  an  analysis  was  presented  which  permitted  the  pre¬ 
diction  of  the  change  of  permeability  as  a  function  of  both  time  and 
position  in  the  porous  medium  while  the  dimensions  of  the  pores  were  being 
changed  by  the  deposition  of  carbon  due  to  thermal  decomposition  of  methane. 

In  the  present  investigation,  the  previous  work  has  been  extended  to  in¬ 
clude  carbon  deposition  and  removal  using  acetylene,  hydrogen,  carbon 
monoxide,  water,  and  methane  for  the  test  gases.  These  were  used  singly 
and  in  mixtures.  The  effects  of  actual  pyrolysis  gases,  obtained  from  a 
phenolic  resin,  were  also  investigated.  The  chars  used  were  produced  from 
cylinders  having  the  BSD  designation  of  R6300-HP. 

Carbon  monoxide  and  hydrogen  were  found  to  have  a  negligible  effect  on  the 
permeability  of  'che  chars.  Both  methane  and  acetylene  decreased  the 
permeability  while  water  increased  the  permeability  of  the  chars.  The 
phenolic  ablation  products  decreased  the  permeability.  It  was  also  noted 
that  the  chars  could  develop  cracks  from  thermal  stresses  alone. 


Tlie  use  of  Che  results  obtained  for  predicting  char  phenomena  during 
actual  re-entry  was  considered.  An  analytical  procedure  was  developed  for 
the  case  of  a  growing  char  in  which  carbon  deposition  is  occurring  due 
to  the  cracking  of  methane. 


A 
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II 


BACKGROUND 


2.1  DARCY’S  LAW 

Efforts  to  correlate  the  flow  rate  of  a  fluid  through  a  porous  material 
with  the  pressure  drop  across  it  have  been  made  for  more  than  a  hundred 
years „ 


The  experiments  of  Darcy  (Reference  2)  and  later  workers  led  to  the  de¬ 
velopment  of  a  differential  equation  which  has  come  to  be  known  aa  "Darcy’s 
law".  For  one  dimensional  flow  this  is 


-  —  ■  H-  y 
dx  K  x 


(2.1) 


vfoere 


P  *  pressure 

x  “  distance  in  x  direction 
K  “  permeability 
p  -  viscosity 

*  gas  velocity  in  x  direction 

Thus,  Darcy 's  law  states  that  the  flow  rate  through  a  porous  medium  is 
proportional  to  the  pressure  graduent  across  the  medium. 

2.2  ANALYTICAL  MODELS  FOR  PERMEABILITY 


In  order  ,to  describe  the  change  in  permeability  with  carbon  deposition  a 
mathematical  model  will  be  needed.  There  has  been  considerable  work  done 
in  predicting  permeability  as  a  function  of  various  structural  parameters 
of  porous  media.  Scheidegger  (Reference  3)  presents  an  excellent  review  of 
the  models  which  have  been  used  to  represent  porous  structures. 


The  simplest  model  of  a  porous  medium  is  that  consisting  of  a  bundle  of 
straight  capillaries.  The  flow,  Q,  through  a  capillary  can  be  described  by 
the  Uagen-Poiseuille  Law: 


Q 


-rr  (54  dP 
128  p  dx 


(2.2) 


where  6  is  the  diameter  of  the  pores  in  the  model. 


By  expressing  this  for  the  case  where  there  are  n  capillaries  per  unit  area 
and  combining  with  Darcy's  law,  one  obtains  the  following  equation  for  the 
permeability: 


K  <= 


mr  6^* 
128 


(2.3) 
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This  equation  can  be  put  into  various  other  font.  Letting  the  porosity 
e  be  represented  by 


1  74 

e  *  jr  n*  fi 


(2.4) 


and  the  specific  internal  surface  area  of  the  porous  matrix  (the  ratio  of 
the  area  of  the  capillaries  to  the  enclosing  voluae  of  the  solid  be  given 

by 

S  ■  ni  6  (2 .5) 

one  can  substitute  for  n  or  6  in  terns  of  e  or  S. 

One  of  the  most  widel't  accepted  theories  for  permeaHlity  was  presented  by 
Kozeny  (Reference  4)  :n  1927.  H2  represented  a  porous  medium  by  channels 
of  various  cross-secti -ns ,  but  of  definite  length.  Kozeny  solved  the  Navier- 
Stokes  equations  (under  certain  assumptions)  for  the  flow  In  channels  and 
then  substituted  his  results  into  Darcy's  law  to  obtain  the  Kozeny  equation: 


(2.6) 


where  e  and  S  are  defined  as  before  and  c  is  the  Kozeny  constant. 


The  Kozeny  equation  contains  an  important  parameter  in  the  surface  area. 

For  this  reason  it  has  been  widely  used  in  determining  surface  areas.  For 
surface  area  determinations,  the  value  of  1/5  for  c  is  generally  assumed. 

2.3  KINETICS  OF  CARBON  FORMATION  AND  REMOVAL 

The  most  important  parameters  for  predicting  the  effect  on  permeability  due 
to  carbon  deposition  or  removal  are  the  rates  of  the  reactions  of 
the  gases  flowing  through  the  char.  The  flow  model  employed  and  the  density 
of  the  deposits  are  of  minor  importance  in  the  analysis  compared  to  the 
kinetics  of  the  reactions  which  occur. 


2.3.1  Methane  Decomposition 

The  kinetics  of  the  decomposition  of  methane  to  carbon  has  been  studied 
extensively.  This  reaction  can  be  represented  by 


W  *  °(«  +  2H2(g, 


(2.7) 


It  is  almost  universally  agreed  that  the  reaction  rate  for  carbon  formation 
from  methane  gas  can  be  represented  by  a  first  order  reaction  of  the  form 


r  ••  k  C 
rc  CH.  CH. 
4  4 


(2.8) 


Palmer  and  Cullis  (Reference  5)  have  combined  the  experimental  findings  on 
the  decomposition  of  methane  to  carbon  of  several  authors  in  an  Arrhenius 


f 


r 


* 


plot.  The  rate  constant  derived  from  this  figure  was  used  in  the  previous 
report  (Reference  1).  It  is 


r 

c 


10 


(15.16  - 


2.28  x  10 


(2.9) 


Several  authors  have  noted  that  hydrogen  inhibits  the  decomposition  of 
methane  in  the  temperature  range  of  1025  to  1500°K.  However,  Kevorkian, 
Heath,  and  Boudart  (Reference  6)  in  the  temperature  range  1660  to  1965*K  and 
Skinner  and  Ruehrwein  (Reference  7)  in  tv  temperature  range  1430  to  1790*K 
report  no  effect  of  hydrogen. 

2.3.2  The  Kinetics  of  Carbon  Formation  from  Acetylene 

Carbon  formation  from  acetylene  has  been  widely  studied  because  of  its  ap¬ 
pearance  in  reactions  involving  carbon  formation.  It  has  been  postulated 
by  Porter  (Reference  8)  in  the  acetylene  theory  of  carbon  formation  that  car¬ 
bon  is  formed  directly  from  acetylene. 

At  lower  temperatures  (below  900°K)  polymerization  is  the  principal  reaction 
of  acetylene.  At  higher  temperatures  the  products  are  carbon,  hydrogen, 
diacetylene,  vinylacetylene,  and  sometimes  methane  (Reference  9  and  8).  As 
a  general  rule  the  last  three  are  formed  in  negligible  quantities  at  the 
high  temperatures . 

Because  of  this,  the  reaction  can  be  represented  by 

C2112  ->  2C  f  H2  (2.10) 

for  the  temperature  regions  of  interest  in  this  investigation. 


Most  investigators  agice  that  the  kinetics  of  acetylene  decomposition  can  be 
represented  by  a  second  order  reaction  of  the  form 


-  r 


c2ii2 


\\  Cc2H2 


(2.11) 


However,  there  are  some  exceptions  to  this  agreement.  Kozlov  and  Knorre 
(Reference  9)  found  their  data  to  be  best  represented  by  a  first  order 
homogeneous  reaction  rate  and  Silcocks  (Reference  10)  found  his  data  to  be 
best  represented  by  a  combination  of  a  first  order  heterogeneous  reaction 
with  the  surface  in  his  reaction  vessel  and  a  second  order  homogeneous 
reaction. 


2.3.3  Carbon  from  Carbon  Monoxide 


The  reaction  describing  this  carbon  formation  is  known  as  the  Boudouard 
re ac Li on: 


2  CO  l  C02  +  C 
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« 


(2.12) 


Kassal  (Reference  11)  has  studied  the  equilibrium  of  this  reaction.  He 
writes  the  reaction  in  the  reverse  direction  so  that  the  equilibrium  con¬ 
stant  is  defined  as  _ 

P 

K  -  (2.13) 

eq  j 

co2 

Table  1  lists  the  values  he  obtained  for  the  equilibrium  constant. 

From  the  equilibrium  data  one  can  see  that  the  reaction  favors  carbon  forma¬ 
tion  only  below  100G*K.  Therefore,  one  might  expect  carbon  to  be  formed  at 
temperatures  below  1000°K.  However,  this  is  not  the  case  since  the  reaction 
rate  becomes  small  at  low  temperatures.  Gaydon  and  Wolfhard  (Reference  12) 
report  that  below  1500°K  the  gas  phase  reaction  does  not  occur.  From  the 
equilibrium  data  one  can  see  that  when  the  reaction  does  occur,  the  amount 
of  carbon  formation  is  negligible. 

2. 3. A  Reaction  of  Carbon  with  Hydrogen 

The  reaction  of  carbon  with  hydrogen  is  the  reverse  of  the  reactions  thus 
far  discussed;  i.e.,  carbon  is  removed  from  the  surface  rather  than  being 
deposited. 

Sanders  (Reference  13)  studied  this  reaction  at  temperatures  from  1825  to 
2550*K.  He  found  that  his  experimental  weight  losses  agreed  with  his 
calculations  based  on  thermodynamic  data  for  the  hydrogen-graphite  theoret¬ 
ical  equilibrium  composition,  with  acetylene  the  major  reaction  product. 

At  lower  temperatures  methane  is  the  major  reaction  product. 

Table  2  gives  the  equilibrium  constants  obtained  from  thermodynamic  data. 
From  the  table  it  can  be  seen  that  at  1800°K  only  methane  and  acetylene 
contribute  measurably  to  the  carbon  weight  loss,  while  at  2600^  only 
acetylene  is  significant. 

The  main  point  to  note  about  the  carbon-hydrogen  reaction  is  that  its 
equilibrium  conversion  is  small  for  the  temperature  range  discussed.  It  is 
not  until  temperatures  in  excess  of  3000®K  that  the  reaction  of  hydrogen 
with  carbon  becomes  significant. 

2,3.5  Reaction  of  Carbon  with  Steam 


As  with  many  heterogeneous  reactious,  the  determination  of  the  kinetics  of 
the  carbon-steam  reaction  poses  many  problems.  Processes  can  be  controlled 
by  gas  film  diffusion,  by  pore  diffusion,  by  the  chemical  reaction  rate, 
and  by  combinations  of  these. 


The  rate  of  the  carbon-steam  reaction  is  often  expressed  as 


•aph2° 


r  - 


1  +  Vh2  +  k6  PH20 


(2.1A) 
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According  co  Ergun  and  Mentser  (Reference  14)  the  individual  reactions 
forcing  the  reaction  sequence  are  a  surface  reaction: 


H2°  +  Cf< 


t 

k. 


2  "2  +  Co 


(2.15) 


(where  C.  and  C  refer  to  free  and  unoccupied  active  carbon  sites  respec¬ 
tively),  followed  by  the  gasification  reaction: 


-&■  CO  +  nC, 


(2.16) 


The  carbon  monoxide  can  also  react  to  form  carbon  dioxide  according  to 


CO  +  C 


r*CO„  -r  C. 


(2.17) 


This  can  be  further  coroolicated  by  the  water-gas  shift  reaction: 


CO  +  11,0  <- 


->Cn,  +  d. 


(2.18) 


Using  the  analysis  of  Ergun  and  Mentser  in  their  review  of  the  carbon-steam 
reaction,  the  rate  of  carbon  removal  can  be  expressed  as 


<  v  r ' 

dn  "4  j  'ct 


dt  K,  +  C  /C„  0 

*<  n_  n- 


(2.15) 


cn 


where  -  - —  is  tiie  rate  of  carbon  removal  fiom  the  surface,  atoms  per 
second;  K  is  the  mass  of  the  solid  carbon;  R^  is  the  equilibrium  constant 
for  reaction  (2.15);  and  k^C^  is  the  product  of  the  rate  constant  for 
reaction  (2.16;  tines  the  total  number  of  active  carbon  sites.  This  ex¬ 
pression  is  on  a  weignt  basis,  however,  a  more  applicable  expression  may 
be  obtained  by  writing  this  in  ter :.\r,  of  the  surface  area,  which  is  the 
actual  governing  parameter  for  tne  reaction  rate,  rather  than  the  weight. 

The  surface  area  per  unit  volume  may  oe  expressed  in  terms  of  the  perme¬ 
ability  using  the  Kozeny  equation  to  obtain 


P~ 

,  R ,  V  k’  C’  li¬ 
do  3  cty  5  K 

_ c _ 

dt  K.  +  C..  /C.  A 
4  lUO 


(2.20) 


where  V  is  the  volume  and  k^  is  the  rate  constant  in  terms  of  surface  area. 
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Ergun  and  Mentser  experimentally  studied  this  reaction  in  the  range  1275 
to  1475°K  and  determined  values  for  and  k^C^.  Their  value  of 
should  hold  for  other  types  of  porous  carbonaceous  systems.  The  term  k^  Cct 
will  give  the  temperature  dependence  of  C^t. 

In  analyzing  their  results,  Ergun  and  Mentaer  determined  that  the  water-gas 
shift  reaction  was  at  equilibrium  at  1475°K  in  their  reaction  vessel  and 
reached  equilibrium  by  the  time  the  1^0  was  502  converted  at  lower  tempera¬ 
tures.  This  fact  can  be  used  for  determining  the  product  distribution  in 
the  reaction  zone. 

2.3.6  Mixtures 

The  effect  of  hydrogen  on  the  rate  of  carbon  formation  from  methane  was 
discussed  in  the  section  on  methana.  There  is  no  information  on  the  effect 
of  hydrogen  on  the  reactions  of  the  other  gases  except  for  the  reaction  of 
water  vapor  with  carbon.  In  this  case  the  effect  of  hydrogen  and  of  carbon 
monoxide  can  be  taken  into  account  with  the  water-gas  shift  equilibrium. 

The  only  other  interactions  that  might  be  expected  are  the  interactions 
between  the  water  vapor  and  the  methane  and  acetylene.  However,  a  litera¬ 
ture  search  has  uncovered  no  information  on  this  subject. 

2.4  DENSITY  OF  DEPOSITED  CARBON 

The  density  of  the  deposit  is  needed  to  determine  the  volume  the  deposited 
carbon  will  occupy  in  the  pores.  The  density  of  deposited  carbon  formed 
by  the  decomposition  of  hydrocarbon  gases  has  been  studied  by  Brown  and 
Watt  (Reference  15).  They  found  that  the  density  of  the  deposited  material 
is  dependent  on  the  temperature  of  deposition.  The  density  varied  from  a 
minimum  of  1.14  g/cro^  at  1975°K  to  a  nuximum  of  2.22  g/cm^  at  2375 °K.  At 
temperatures  from  2300  to  2400°K  and  higher,  the  density  was  essentially 
constant.  Figure  1  shows  the  change  in  density  with  temperature  of 
deposition. 


-8- 


Ill 


EQUIPMENT 


3.1  INDUCTION  FURNACE 

The  induction  furnace  for  heating  the  samples  to  temperatures  character¬ 
istic  of  ablation  is  discussed  in  detail  in  Reference  1.  The  apparatus 
was  capable  of  operating  at  temperatures  up  to  3000*K  for  indefinite 
periods  of  time.  It  was  leakproof  so  that  air  could  not  enter  the  system 
and  the  gases  being  passed  through  the  system  could  not  escape.  Provision 
was  also  made  for  accurate  optical  pyrometer  temperature  measurements  de¬ 
spite  deposition  of  opaque  matter  on  the  quartz  containing  walls. 

The  furnace  without  the  radiation  shield  is  shown  in  Figure  2.  As  indi¬ 
cated  by  the  name,  induction  heating  was  used  to  achieve  the  necessary  tem¬ 
perature  levels.  The  power  was  supplied  by  a  Lepel  Model  T25-3,  25  kilo¬ 
watt  generator  with  a  frequency  range  from  180  kc  to  450  kc. 

The  gases  used  in  the  experiments  were  introduced  into  the  specimen  cylinders 
through  an  inlet  in  the  bottom  of  the  furnace.  A  thin  pyrolytic  graphite 
tube  was  used  to  carry  the  gases  from  the  bottom  of  the  furnace  to  the  hole 
in  the  bottom  sample  holders.  The  gases  then  flowed  radially  outward 
through  the  specimen  and  exited  through  the  top  of  the  furnace. 

A  Leeds  and  Northrup  8630-series  optical  pyrometer  was  used  to  measure  the 
surface  temperatures  of  the  chars  during  the  transpiration  experiments. 

This  was  a  single-adjustment,  potentiometer- type  of  instrument.  The  sample 
was  viewed  through  a  quartz  sight  tube  (1/4  inch  O.D.)  which  was  attached 
to  the  inside  surface  of  the  quartz  cylinder  about  1  1/2  Inches  below  the 
level  of  the  specimen.  Part  of  the  cooling  gas  was  diverted  to  flow  past 
the  sight  tube  to  keep  the  sight  path  clear.  A  mirror  was  placed  on  the 
bottom  steel  end  plate  to  reflect  the  radiation  coming  through  the  sight 
tube  at  an  angle  more  convenient  for  viewing  with  the  pyrometer. 

In  order  to  obtain  the  high  temperatures  desired  in  the  furnace,  it  was 
necessary  to  use  a  radiation  shield.  A  pyrolytic  graphite  cylinder  (1  1/2 
inches  high  and  approximately  1/4  inch  thick)  was  used  for  this  purpose. 

The  graphite  was  oriented  so  that  the  direction  of  lowest  conductivity  was 
in  the  radial  direction. 

3.2  FURNACE  MODIFICATION  FOR  WATER  VAPOR 

The  gas  inlet  system  to  the  furnace  had  to  be  modified  to  incorporate  water 
vapor  as  the  test  gas.  This  was  accomplished  by  injecting  the  water  through 
a  small  nozzle  into  the  primary  test  gas  stream  where  the  water  was  vapor¬ 
ized.  The  nozzle  was  made  from  a  stainless  steel  Luer  syringe  27  gauge 
needle.  The  primary  gas  stream  was  heated  to  approximately  120*C  by 


-9- 


wrapping  the  gas  line  with  electrical  heating  tape.  A  Variac  was  used  to 
control  the  temperature  of  the  tape.  The  liquid  water  flow  was  monitored 
with  a  rotameter.  The  water  injection  system  is  depicted  in  Figure  3. 

3.3  FURNACE  MODIFICATION  FOR  YARN 

The  yarn  consisted  of  a  phenolic  resin  impregnated  in  a  carbon  yarn  which 
was  then  cured.  The  yarn  was  passed  through  the  center  of  the  char  in  the 
furnace  where  the  phenolic  resin  was  partially  pyrolyzed.  The  system  used 
for  handling  the  yam  is  shown  in  Figure  4.  Prior  to  an  experiment,  the 
yarn  was  wound  on  the  source  spool  and  then  threaded  through  the  gas  inlet 
to  the  furnace  and  back  down  to  a  motor  driven  wind  up  spool  where  the  yarn 
was  rewound  after  passing  through  the  furnace.  The  quartz  tube  on  the  top 
char  holder  was  necessary  to  provide  a  large  smooth  turning  radius  for 
the  yam,  since  the  yam  could  not  sustain  sharp  turns.  The  motor  and 
spools  were  contained  in  a  pressurized  tank;  the  test  gas  passed  through  the 
tank  and  into  the  furnace.  The  motor  speed  was  controlled  with  a  Variac. 

3.4  PORE  SIZE  DISTRIBUTION 

A  1960  model  Amin co-Wins  low  Porosimeter  was  used  to  determine  the  pore  size 
distribution.  This  was  a  5000  psi  hand-operated  model  which  could  be  used 
to  determine  the  volume  of  pores  ranging  in  size  from  100  to  0.035  microns. 
The  results  were  semiquantitative,  since  the  diameter  of  any  pore  as  de¬ 
termined  by  the  porosimeter  was  the  smallest  diameter  through  which  the 
mercury  must  pass  to  enter  the  pore.  This  is  known  as  the  "bottleneck 
effect".  The  data  obtained  from  the  mercury  porosimeter  were  the  total 
volume  of  pores  that  are  filled  with  mercury  at  given  pressures.  From  this 
information  one  could  calculate  a  capillary  pressure  curve  which  in  turn 
could  be  related  to  the  pore  size  distribution. 

3 .5  POROSITY 

The  porosity  was  determined  by  a  Boyles'  law  porosimeter  which  uses  helium 
gas.  In  this  apparatus  a  metered  amount  of  helium  was  brought  into  the 
sample  chamber.  By  knowing  the  total  volume  of  the  sample  chamber  and 
measuring  the  pressure  before  and  after  the  gas  was  in  the  apparatus,  one 
could  determine  the  volume  of  the  solids  and  closed  pores,  and  hence  the 
porosity.  The  pressure  was  measured  by  a  fused  quartz  precision  pressure 
gauge  manufactured  by  the  Texas  Instrument  Company. 

3.6  ANALYTICAL  BALANCE 

All  sample  weighings  were  performed  on  an  automatic  balance  manufactured 
by  Wm.  Ainsworth  &  Sons,  Inc.  The  balance  is  capable  of  weighing  to 
0.001  grams.  These  weighings,  in  conjunction  with  the  helium  porosimeter 
results,  were  used  to  determine  real  and  apparent  densities. 
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3-7  PRESSURE  :;EASURr;:;^.f 

The  pressure  drop  across  the  sample  was  also  monitored  with  the  quartz 
pressure  gauge.  The  gauge  output  was  in  digital  form;  a  strip  recorder 
was  used  to  obtain  a  rerord  of  the  pressure  variation  in  the  course  of 
an  experiment. 

3.8  FLOW  RATES 

Rotameters  were  used  for  monitoring  the  flaw  rates  of  the  gases.  These 
were  Brooks  rotameters  and  had  maximum  flow  rates  ranging  from  51  cc/min. 
to  24.4  SCr‘I  of  air. 

3.9  GAS  ANALYSIS 

The  effluent  gases  from  the  furnace  were  analyzed  with  a  gas  chromato¬ 
graph.  The  analysis  did  not  give  the  exact  composition  of  the  gases 
leaving  the  sample  since  the  gases  continued  to  react  after  that  time. 
This  effect  was  minimized  by  a  side  effect  of  the  gas  used  to  keep  the 
sight  tube  clear.  This  gas  tended  to  quench  the  reaction  by  cooling  the 
effluent  gases  immediately  upon  their  exit  from  the  sample.  A  Beckman 
GC-5  chromatograph  equipped  with  a  temperature  programmer  was  used  during 
most  of  the  experimental  work. 
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IV 


CHARACTERI ZATICH  OF  CHARRED  CYLINDERS 


4.1  DESCRIPTION  OF  CYLINDERS  AS  SUPPLIED 


Eight  charred  cylinders  were  supplied  by  Avco/RAD.  All  cylinders  were  made 
from  compression  molded  high  purity  RAD-6300-HP  carbon  phenolic  material 
using  the  dixie  cup  method  with  a  30  degree  layup  angle  to  the  pipe  axis. 
Before  charring,  the  cylinders  were  twelve  inches  in  length  with  an  outside 
diameter  of  1.25  inches  and  a  wall  thickness  of  0.25  inches. 


Btu 


Charring  was  accomplished  by  passing  plasma  gases  from  the  Avco/RAD  Ten 
Megawatt  Arc  Facility  through  the  inside  of  the  cylinders  for  periods 
ranging  from  14  to  27  seconds.  The  heat  transfer  rate  to  the  cylinders  was 
about 

1100 

ft  sec. 

The  peak  backface  temperatures  ranged  from  1120  to  2000°F.  The  thickness 
measured  at  the  center  of  the  cylinders  varied  from  0.08  to  0.19  inches 
after  charring.  The  test  conditions  for  the  charring  are  summarized  in 
Table  3. 


It  should  be  emphasized  that  considerable  variation  of  structural  properties 
were  evident  in  the  angular  direction  at  any  given  axial  position  in  the 
cylinders  as  supplied.  This  was  apparently  due  to  the  non-symmetry  of  the 
plasma  jet  used  in  the  charring  process. 


4.2  CHARACTERIZATION  OF  THE  CYLINDERS  AS  SUPPLIED 

The  cylinders  were  cut  Into  one  inch  long  segments.  The  location  of  these 
segments  from  the  hot  end  during  the  charring  process  Is  designated  in  the 
sample  name.  The  dimensions  of  the  samples  used  are  listed  in  Table  4. 

Also  shown  are  the  real  and  apparent  densities,  the  porosity,  the  perme¬ 
ability,  and  the  average  pore  diameters. 

Figure  5  is  a  photomicrograph  cf  an  as-received  char.  Significant  deposi¬ 
tion  may  be  noted  near  the  inside  surface.  Almost  no  deposition  is  evident 
near  the  outside  surface, 

4.3  COMPARISON  WITH  COMMERCIAL  CARBONS  AND  GRAPHITES 

As-received  characterization  of  commercial  carbons  and  graphites  is  shown  in 
Table  5  for  comparison  purposes.  The  carbons  and  graphites  were  obtained 
from  Union  Carbide  Corporation. 
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V 


CONDITIONING  OF  THE  CHARS 


The  chars  in  the  as-received  condition  contained  unablated  residue.  This 
was  eliminated  by  passing  an  inert  gas  through  the  char  at  high  tempera¬ 
tures  in  the  furnace.  The  inert  gas  used  was  helium.  Al1  chars  used  in 
deposition  runs  were  conditioned  in  this  manner.  The  characterization 
results  after  conditioning  are  summarized  in  Table  6.  It  can  be  seen  that 
the  real  density  after  conditioning  was  about  4%  higher  than  before  con¬ 
ditioning.  The  porosity  of  the  samples  was  typically  28%  before  condition¬ 
ing  and  32%  after  conditioning.  In  general,  the  after  conditioning  per¬ 
meabilities  were  about  one  and  half  times  larger  than  the  before  condition¬ 
ing  permeabilities.  Most  of  the  chars  were  conditioned  at  2600°K  for 
40  minutes.  The  time  for  conditioning  at  the  highest  temperature  and  the 
highest  temperature  are  both  listed  in  Table  6. 

It  is  of  interest  to  compare  the  permeabilities  of  the  high  purity  chars 
used  in  the  present  study  with  the  chars  of  Reference  1  and  also  with  the 
commercial  carbons  and  graphites.  Figure  6  shows  such  a  comparison.  It 
is  immediately  obvious  that  the  permeabilities  of  the  present  chars  are 
an  order  of  magnitude  lower  than  the  previous  chars  and  several  orders  of 
magnitude  lower  than  the  carbons  and  graphites.  All  of  the  chars  shown  in 
Figure  6  were  conditioned  at  about  2600°K  for  33  minutes  or  more. 

The  pore  size  distribution  of  a  char  before  conditioning  is  compared  with 
a  char  after  conditioning  in  Figure  7.  Conditioning  had  only  a  small  effect 
on  the  pore  size. 

Three  chars  were  conditioned  at  high  temperature,  about  2900oK.  All  three 
chars  cracked  within  ten  minutes  at  this  temperature.  None  of  these  chars 
was  cemented  in  the  sample  holders  for  conditioning.  After  two  chars  had 
cracked  at  high  temperature,  great  care  was  taken  with  the  third  sample 
(R6300-HP-DC4-0 .6FH-XS)  to  determine  the  reasons  for  cracking.  The  sample 
holders  were  remachined  to  obtain  clean  smooth  surfaces.  The  top  (spring 
loaded)  support  was  not  used  in  order  to  minimize  unnecessary  mechanical 
forces  on  the  char.  The  top  holder  was  needed  to  prevent  radiation  losses; 
it  weighed  less  than  25  grams. 

The  char  was  conditioned  for  ten  minutes  at  26fj0°K  without  difficuJty. 

The  temperature  was  then  increased  to  2875°K.  After  less  than  4  minutes, 
the  temperature  in  the  char  dropped  rapidly  indicating  a  crack  in  the  char. 
When  it  was  removed  from  the  furnace,  the  crack  was  confirmed;  a  picture 
of  the  crack  is  shown  in  Figure  8.  A  large  pit  on  the  inside  surface  also 
developed  during  conditioning  and  may  be  seen  in  Figure  9.  The  depth  of 
this  pit  was  about  .07".  Since  there  were  no  mechanical  forces  on  the 
sample,  the  crack  must  be  ascribed  to  thermal  stresses. 
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A  photomicrograph  of  a  char  after  conditioning  at  2660°K  is  shown  in 
Figure  10.  Comparing  this  with  a  sample  as  received,  the  after-condition- 
ing  sample  is  more  uniform  in  the  radial  direction.  Also  there  are  dark 
deposits  visible  at  the  outside  surface  similar  to  those  at  the  inside 
surface.  These  deposits  at  Che  outside  surface  are  a  result  of  the  con¬ 
ditioning  process.  They  are  deposited  carbon,  formed  during  conditioning 
by  the  decomposition  of  the  gases  coming  from  the  unablated  material. 

The  flow  direction  of  the  gas  was  from  the  inside  out,  so  that  the  deposi¬ 
tion  occurred  in  the  outer  section  of  the  sample. 
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VI 


CARBON  DEPOSITION  AND  REMOVAL 


6.1  EXPERIMENTAL  PROCEDURE 


To  simulate  the  char  densification  that  occurs  in  the  char,  gases  typical 
of  the  pyrolysis  products  of  phenolic  resins  CO,  i^O,  and  CH^) 

both  singly  and  in  mixtures  were  passed  through  cnar  samples  while  the 
samples  were  inductively  heated.  The  effects  of  these  gases  on  the  char 
were  then  determined. 


For  test  gases  that  could  be  obtained  in  compressed  cylinders  (He,  CH^, 

C2H2,  and  CO)  the  general  experimental  procedure  was  to  heat  the  sample 
to  the  desired  temperature  using  helium  only.  The  test  gas  was  then  allowed 
to  flow  through  the  sample.  Pressures  and  flow  rates  were  monitored,  and 
from  these  permeability  as  a  function  of  time  wai  calculated. 


-> 

The  procedure  for  water  vapor  was  the  same  as  that  above,  except  that  the 
water  stayed  in  the  liquid  form  until  after  it  was  mixed  with  the  other 
gases .  It  was  then  vaporized  and  entered  the  furnace  as  a  gas . 


n 

For  experiments  with  yarn,  the  yam  was  first  dried  in  a  vacuum  oven  to  re¬ 
move  adsorbed  water  vapor.  It  was  then  placed  in  the  experimental  appa¬ 
ratus.  The  motor  pulling  the  yam  through  the  furnace  was  started  and  the 
temperature  was  immediately  raised  to  the  teat  temperature.  The  weight 
loss  of  the  yam  was  jdetermined  by  weighing  the  yam  that  had  passed  through 
the  furnace  and  also  a  control  piece  that  was  dried  but  had  not  passed 
through  the  furnace.  By  measuring  the  lengths  of  the  two  pieces,  the  weight 
loss  per  unit  length  could  be  determined,  and  from  this  the  total  weight 
loss . 


6.2  CALCULATION  ,0F  PERMEABILITIES  FROM  THE  EXPERIMENTAL  DATA 


When  the  experimental  data  were  reduced  to  obtain  permeabilities,  the  change 
in  gas  volume  due  to  reaction  was  neglected  so  as  to  avoid  undue  complex¬ 
ities  in  the  computations .  For  the  relatively  small  concentrations  of  re¬ 
acting  gases  .used  in  the  experiments  this  assumption  would  not  cause  any 
appreciable  error. 

In  determining  the  mobility  (the  ratio  of  permeability  to  initial  perme¬ 
ability),  the  initial  value  of  the  permeability  used  was  the  value  obtained 
8 lightly  after  the  test  gas  was  started  extrapolated  to  time  zero.  Thus 
errors  caused  by  neglecting  the  reaction  of  the  test  gas  cancelled  each 
other  in  calculating  mobility,  and  any  changes  in  permeability  that  occurred 
during  startup  operations  were  not  erroneously  included  in  the  results. 

The  permeabilities  after  deposition  in  the  chars  are  summarized  in  Table  7; 
the  carbon  and  graphite  results  are  shown  in  Table  8.  A  comparison  of  the 
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apparent  density,  permeability,  and  porosity  of  the  chars  in  the  as- 
received,  after  conditioning  and  after  deposition  conditions  may  be 
found  in  Table  9,  During  the  experiments  the  temperature  and  flow  rates 
varied  somewhat.  In  presenting  the  results  in  the  tables  and  figures, 
average  values  were  used. 

6.3  RESULTS  USING  METHANE 

The  analysis  presented  in  the  previous  report  for  the  prediction  of  the 
change  in  permeability  as  a  function  of  time  and  position  for  carbon  deposi¬ 
tion  by  the  thermal  decomposition  of  methane  has  been  verified. 

Complete  data  and  results  for  deposition  in  carbon  and  graphite  samples  can 
be  found  in  the  dissertation  by  Schvind  (Reference  16) . 

A  comparison  of  the  analysis  with  the  experimental  results  from  two  of  the 
runs  is  shown  in  Figures  11  and  12.  The  results  are  plotted  in  terms  of 
mobility.  It  is  seen  that  the  analysis  provides  an  accurate  tomans  of  pre¬ 
dicting  the  effect  of  carbon  deposition  from  the  decomposition  of  methane 
in  these  types  of  porous  media. 

Figure  13  shows  the  effect  of  methane  deposition  on  the  permeability  (mo¬ 
bility)  of  a  char  sample.  Some  of  the  methane  decotqioaed  prematurely  during 
this  experiment.  This  can  be  seen  in  Figure  14  which  is  a  photomicrograph 
of  the  sample  after  the  run.  The  deposit  at  the  inside  surface  is  the 
material  just  slightly  ciifferent  in  color  from  the  epoxy  resin  that  was 
used  for  impregnating  the  sample  in  preparing  it  for  photomicrography.  The 
premature  decomposition  was  due  to  the  low  permeability  of  the  char  and  the 
resulting  low  flow  rate. 

6.4  RESULTS  USING  ACETYLENE 

Acetylene  decreases  the  permeability  of  the  chars  by  the  deposition  of 
carbon.  The  decrease  in  permeability  is  faster  than  that  observed  with 
methane.  Some  of  the  experimental  results  are  presented  in  Figures  15 
through  19. 

Figures  20  and  21  are  photomicrographs  of  the  R6300-HP-DC1-8.2FH-XS  sample 
after  carbon  deposition  using  acetylene.  Here  again,  some  of  the  acetylene 
decomposed  prematurely.  This  can  be  seen  as  a  deposit  at  the  inside  surface 
with  little  deposition  in  the  pores  of  the  sample  in  Figure  21.  Figure  20 
is  a  different  view  of  the  same  sample.  In  this  region  of  the  sample  the 
acetylene  did  not  decompose  prematurely  and  deposit  at  the  inside  surface 
but  in  the  sample,  the  black  spots  in  Figure  20  being  deposited  carbon. 

6.5  RESULTS  USING  CARBON  MONOXIDE 

In  Section  II  it  was  predicted  on  Che  basis  of  equilibrium  data  that  when 
carbon  monoxide  reacts  to  form  carbon,  the  amount  of  carbon  formed  will  be 
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negligible.  One  sample,  R6300-HP-DC1-9 .4FH-XS,  was  tested  with  carbon 
monoxide  at  three  temperatures  (1620°,  2040®,  and  2760®K) .  At  all  three 
temperatures  the  permeability  remained  constant;  the  results  agreed  with 
the  prediction. 

6.6  RESULTS  USING  HYDROGEN 

In  Section  II  it  was  also  predicted  that  the  hydrogen  reaction  with  carbon 
is  equilibrium-limited.  Not  until  temperatures  are  in  excess  of  3000°K 
would  it  be  expected  that  the  effect  of  hydrogen  becomes  important  in  chars. 
The  R6300-HP-DC1-4 .2FH-XS  sample  was  used  for  runs  with  hydrogen  at  1490*. 
1995°,  2660®K.  The  R6300-HP-DC1-5 . 4FH-XS  sample  was  used  for  runs  with 
hydrogen  at  1505°  and  2400°K.  Except  at  the  two  high  temperature  values 
there  was  no  effect  on  the  char  permeability. 

At  the  temperature!  of  2400°  and  266D°K.  both  samples  cracked.  It  is  sus¬ 
pected  that  the  cracking  occurred  because  of  thermal  stresses  for  the  3ame 
reason  that  cracking  occurred  during  conditioning. 

For  the  samples  that  didn't  crack,  the  experinKnts  agree  with  what  one  would 
predict  using  equilibrium  data.  That  is,  hydrogen  has  a  negligible  effect 
on  the  permeability  for  temperatures  up  to  2000CK  and  probably  for  tempera¬ 
tures  up  to  3000°K. 

6.7  RESULTS  USING  WATER  VAPOR 

The  effect  of  water  on  the  chars  is  to  Increase  their  permeability.  The 
water  reacts  with  the  carbon  surfaces  in  the  pores  tc  remove  carbon  from 
the  char. 

Figure  22  shows  the  effect  of  water  on  carbon  chars.  The  permeability  in¬ 
creases  with  time.  Figure  23  is  a  photomicrograph  of  one  of  the  samples 
after  exposure  to  water.  Comparing  this  with  a  sample  after  conditioning. 
Figure  10,  one  can  see  the  dramatic  erosicn  effect  of  the  water.  The  in¬ 
side  section  of  the  sample  has  little  solid  material  remaining.  This  ef¬ 
fect  of  the  water  can  also  be  seen  in  the  pore  size  distribution,  Figure  24. 

6.8  RESULTS  USING  MIXTURES 

Several  authors  have  noted  an  inhibition  of  the  rate  of  methane  decomposi¬ 
tion  by  hydrogen  at  relatively  low  temperatures.  Me thane -hydro gen  mixtures 
were  used  in  runs  with  two  samples.  The  results  are  shewn  in  Figures  25 
and  26.  To  determine  the  effect  of  hydrogen,  the  two  experiments  were 
compared  with  the  analysis  for  methane  alone,  assuming  that  the  hydrogen  was 
an  inert  gas.  The  two  experiments  using  methane  are  compared  with  Schwind's 
experimental  data  in  Figure  27.  The  data  are  plotted  in  terms  of  the  log 
of  the  experimental  rate  constant  for  methane  decomposition  versus  recipro¬ 
cal  temperature.  From  this  figure  it  can  be  seen  that  for  temperatures 
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greater  than  1580°K,  hydiogen  does  not  affect  the  decomposition  of  methane. 

The  deviation  of  Schwind's  experimental  data  from  his  theoretical  line  at 
low  temperature  was  due  to  the  fact  that  at  the  lew  temperatures  the  carbon 
deposition  was  a  result  of  the  decomposition  of  the  impurities  in  his 
methane  gas  rather  than  a  result  of  the  decomposition  of  the  methane.  For 
temperatures  greater  than  2100°K  the  methane  decomposed  before  actually 
entering  the  sample.  Since  the  same  purity  of  methane  gas  was  used  for  the 
two  methane-hydrogen  experiments  that  Schwind  used,  the  line  that  best  fit 
his  data  was  used  to  test  the  effect  of  the  methane-hydrogen  mixture. 

Figures  28  and  29  show  the  results  using  acetylene  and  hydrogen.  With  the 
char  sample  the  acetylene  decomposed  prematurely.  However,  using  the  re¬ 
sults  of  the  acetylene  and  hydrogen  with  a  carbon  sample  and  comparing  them 
with  results  with  just  acetylene,  it  can  be  seen  that  hydrogen  does  not 
affect  the  carbon  deposition  from  acetylene.  Figures  30  and  31  show  the 
results  for  acetylene-methane  and  acetylene -carbon  monoxide  mixtures. 

Figures  32  through  35  show  the  effect  of  carbon  deposition  and  removal 
occurring  simultaneously.  By  varying  the  concentrations  of  the  gases,  the 
permeability  can  be  made  to  increase  or  decrease.  By  increasing  the  concen¬ 
tration  of  the  gases  that  deposit  carbcn,  CII^  and  C2H2,  the  permeability 
can  be  made  to  decrease  more  rapidly.  By  increasing  the  concentration  of 
the  gas  which  removes  carbon  from  the  pores,  ^0,  the  permeability  can  be 
made  to  increase  more  rapidly. 

6 .9  RESULTS  USING  YARNS  IMPREGNATED  WITH  PHENOLIC 

Three  experiments  were  made  using  yarns  impregnated  with  phenolic  so  that 
the  effect  of  actual  products  of  ablation  on  the  permeability  of  chars  could 
be  studied.  In  all  cases  the  effect  of  the  decomposition  products  of  the 
yarn  was  to  decrease  the  permeability  of  the  chars.  The  yam  used  was  Union 
Carbide  Co.  VYB  70  1/2,  impregnated  with  Monsanto's  resin  SC  1008.  It  was 
prepared  by  Lockheed  Missile  and  Space  Company.  The  resin  content  of  the 
yam  was  23.9  weight  per  cent. 

The  yarn  was  quite  brittle  and  required  gentle  handling.  In  the  design  of 
the  experimental  apparatus  care  was  taken  to  insure  that  the  yam  did  not 
have  to  negotiate  any  sharp  turns .  The  maximum  temperature  to  which  the 
yam  could  be  exposed  and  still  be  handled  satisfactorily  was  1500°K. 

Before  the  yam  was  used  in  the  furnace  it  was  treated  in  a  vacuum  oven  to 
remove  adsorbed  water  vapor.  The  adsorbed  water  vapor  amounted  to  more 
than  10%  of  the  total  weight  of  the  yam. 

The  R6300-HP-DC8-5 .0FH-XS  char  sample  was  tested  with  the  yam  in  conjunc¬ 
tion  with  a  helium  flow  rate  of  0.108  ft3/min  at  1265°K  using  56.3  feet  of 
yam.  The  mobility  decreased  from  one  to  0.642  in  4.3  minutes. 

The  R6300-HP-DC8-0.6FH-XS  sample  was  tested  with  the  yam  with  a  helium  flow 
rate  of  0.148  ft^/min  at  1265*K  using  91.7  feet  of  yam.  During  this 
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experiment,  leaks  developed  in  the  cement  holding  the  sample.  Even  with 
the  leaks  the  mobility  decreased  to  0.695  in  6.3  minutes. 

The  R6300-HP-DC8-3.9FH-XS  sample  was  tested  with  yarn  with  a  helium  flow 
rate  of  0.120  ft3/min  at  1440°K  using  97.4  feet  of  yam.  The  mobility 
decreased  to  0.578  in  6.6  minutes. 

Exact  determination  of  the  resin  weight  loss  of  the  yam  could  not  be  made. 
Due  to  the  brittleness  of  the  yam,  the  yam  could  break  at  any  time  during 
an  experiment  so  that  the  length  of  yam  to  be  used  in  the  furnace  was  not 
known  beforehand.  This  meant  that  average  values  of  the  weight  per  unit 
length  had  to  be  used  to  determine  the  weight  loss.  At  the  temperatures 
used  the  weight  loss  of  the  resin  was  small  and  of  the  same  order  of  magni¬ 
tude  as  the  variation  in  the  weight  of  the  yarn.  After  the  yam  was  charred, 
it  was  cut  into  15  foot  lengths,  several  of  which  had  weights  greater  than 
an  average  value  for  an  uncharred  piece. 

Another  problem  in  determining  the  weight  loss  was  the  large  amount  of 
water  that  could  be  adsorbed  by  the  sample.  During  the  time  the  yam  was 
exposed  to  the  air  after  drying  (while  being  transferred  from  the  vacuum 
oven  to  the  induction  furnace  and  while  being  transferred  from  the  induc¬ 
tion  furnace  to  the  gravimetric  balance)  some  water  could  be  adsorbed  by 
the  yam.  For  all  the  experiments  the  resin  weight  loss  was  les3  than 
0.005  gm/ft. 

Only  the  final  mobility  for  the  yarn  experiments  is  presented  rather  than 
a  plot  of  mobility  as  u  function  of  time.  This  was  done  due  to  the  rel  i- 
tive  values  of  the  volume  of  the  yam  tank  as  compared  to  the  flow  rate  of 
the  helium.  The  actual  pressure  lagged  behind  the  true  pressure  due  to  the 
time  required  to  change  the  pressure  in  the  tank  using  the  low  flow  rate. 
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VII 


CHAR  GROWTH  SIMULATION 


As  vai  described  in  the  Introduction ,  the  ultimate  aim  of  the  work  which 
was  performed  under  this  and  the  previous  contract  was  to  provide  a  solid 
basis  for  calculations  of  the  change  of  permeability  and  pressure  drop 
across  a  char  growing  under  re-entry  conditions .  It  is  believed  that  this 
aim  has  been  accomplished,  even  though  the  development  of  the  detailed 
calculational  routines  was  not  part  of  the  planned  work  effort.  A  simpli¬ 
fied  model  of  a  growing  char  was  constructed  for  the  case  where  a  methane- 
inert  mixture  is  the  gas  flowing  through  the  pores  in  the  char,  from  the 
reaction  zone  to  the  outside  surface.  The  equations  describing  carbon  de¬ 
position  for  this  situation,  which  had  been  developed  previously  (Refer¬ 
ence  l),  were  used  for  this  work.  The  calculations  were  intended  to 
demonstrate  the  feasibility  of  using  the  type  of  data  obtained  in  the  course 
of  the  project  for  predicting  phenomena  more  closely  associated  with  actual 
re-entry  situations. 

The  model  employed  is  a  simple  model  designed  to  show  how  char  densifica- 
tion  affects  the  permeability  of  the  char.  For  this  reason  the  rate  of 
char  growth  and  temperature  profile  were  arbitrarily  set  with  the  emphasis 
of  the  work  being  on  the  effect  of  the  carbon  formation. 

For  the  present  example  a  linear  char  growth  rate  was  assumed  with  the  char 
reaching  a  thickness  of  .05  inches  in  10  seconds ,  A  linear  temperature 
profile  was  also  used  with  the  temperature  at  the  reaction  plane  being 
set  at  900*K  and  the  temperature  at  the  outside  surface  for  the  maximum 
char  thickness  being  2500eK.  It  was  assumed  that  the  phenolic  decomposed 
at  a  reaction  plane  rather  than  throughout  a  zone  of  reaction,  with  the 
products  of  the  decomposition  being  methane  and  argon.  The  initial  char 
permeability  was  set  at  10“**  ft2,  the  mass  flux  of  methane  at  .469  x  10"2 
lbm/ft2  tec  and  the  mass  flux  of  argon  at  7.04  x  1Q“2  lbm/ft2  sec.  These 
conditions  were  not  meant  to  be  necessarily  representative  of  the  condi¬ 
tions  in  an  actual  re-entry  situation.  They  were  used  only  to  illustrate 
the  type  of  simulation  that  could  be  performed  if  the  actual  re-entry  con¬ 
ditions  were  used. 

Figure  37  presents  the  effects  of  carbon  deposition  occurring  in  the  char 
pores  in  a  simulated  re-entry  condition.  As  in  an  actual  re-entry  situa¬ 
tion  the  char  thickness  increases  with  time.  The  gaseous  products  (in 
this  case  methane  and  argon)  coming  from  the  pyrolysis  of  the  resin,  flow 
through  the  char,  decomposing  to  form  carbon  in  the  pores  of  the  char. 

It  is  seen  from  Figure  37  that  the  carbon  deposition  does  not  have  an  ap¬ 
preciable  effect  on  the  mobility  of  the  outside  surface  of  the  char  until 
after  six  seconds  when  the  surface  temperature  reaches  1850°K.  It  is  not 
until  this  temperature  ia  reached  that  the  rate  of  decomposition  of  the 
methane  becomes  significant  compared  to  the  residence  time  of  the  methane 
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in  the  char.  The  mobility  versus  time  curve  becomes  constant  because  all 
the  methane  has  decomposed  by  the  time  it  has  reached  the  outside  surface 
for  times  greater  than  9  seconds.  The  pressure  drop  from  the  reaction 
plane  to  any  point  in  the  char,  the  outside  surface  temperature,  and  char 
thickness  are  also  plotted  for  the  assumed  conditions. 

The  program  for  this  solution  is  in  the  Appendix. 

The  general  outline  for  this  program  could  be  utilized  in  developing  a 
more  realistic  model,  including  the  actual  gas  mixtures  involved,  varia¬ 
tions  in  the  temperature  gradient  through  the  char  with  time,  etc.  This 
should  permit  prediction  of  the  properties  of  the  char  (permeability, 
porosity,  and  density)  and  the  composition,  pressure,  and  flow  rate  of  the 
gaseous  ablation  products  as  a  function  of  position  in  the  char  and  of  time. 
This  information,  coupled  with  structural  analyses  of  the  char,  would  make 
it  possible  for  the  designer  to  determine  what  combination  of  conditions 
might  lead  to  mechanical  failure  of  the  char.  He  would  then  be  able  to 
specify  variation  of  such  properties  as  weave  of  the  reinforcing  cloth, 
polymer  composition,  and  shield  thickness  in  order  to  alter  the  performance 
of  the  char,  if  desired. 

It  should  be  noted,  however,  that,  for  the  reasons  already  discussed,  exper¬ 
imental  data  on  chars  alone,  as  obtained  in  the  present  work,  are  not 
sufficiently  accurate  to  serve  as  a  reliable  base  for  this  type  of  model. 

It  is  believed  that  a  combination  of  information  from  more  uniform  (and 
more  permeable)  porous  carbonaceous  materials  and  from  the  chars  is  re¬ 
quired  for  the  analysis.  To  be  more  specific,  results  for  the  deposition 
kinetics  of  gases  and  gas  mixtures  obtained  on  the  more  uniform  (and  re¬ 
producible)  materials  could  be  combined  with  data  on  the  pore  size  distribu¬ 
tion  and  internal  surface  area  of  heat  shield  materials  in  various  stages 
of  charring  to  yield  fairly  accurate  results.  The  term  "fairly  accurate" 
must  be  inserted  here  because  of  the  inherent  inhomogeneities  and  non¬ 
uniformities  of  all  practical  ablative  composites,  as  well  as  the  difficulty 
in  predicting  external  effects . 

A  necessary  input  for  the  exact  char  growth  analysis  would  be  the  heat  flux 
at  the  surface  as  a  function  of  time,  the  composition  of  the  original 
gaseous  ablation  products,  the  rate  of  decomposition  of  the  polymer  as  a 
function  of  temperature,  the  thermal  conductivity  and  density  of  the  char, 
and  the  permeability  or  the  initial  char  material  formed.  With  this  in¬ 
formation  the  rate  of  char  growth,  the  char  thickness,  the  temperature 
profile,  the  gas  flow  rate,  the  porosity,  the  density,  the  pressure  and  the 
permeability  of  the  char  could  be  calculated  as  a  function  of  time  ar.d  po¬ 
sition.  Also  the  assumption  used  in  the  simplified  char  growth  simula¬ 
tion  program  (with  methane),  of  the  polymer  decomposing  at  a  plane  rather 
than  throughout  a  zone,  could  be  eliminated. 
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VIII 


CONCLUSIONS  AND  RrXOM.ll_\'l)Ai  IONS 


The  effect  of  carbon  deposition  and  removal  in  ablative  chars  was  studied, 
using  gases  typical  of  the  ablation  products  of  phenolic  decomposition  both 
singly  and  in  mixtures.  It  was  found  that,  by  themselves,  carbon  monoxide 
and  hydrogen  had  a  negligible  effect  on  the  permeability  of  chars.  Methane 
and  acetylene  decreased  the  permeability,  the  effect  of  acetylene  being 
more  rapid  than  methane.  Hydrogen  did  not  affect  the  rate  of  carbon  forma¬ 
tion  from  either  methane  or  acetylene,  it  acted  only  as  a  diluent.  Water 
increased  the  permeability  of  chars.  The  composition  of  mixtures  of  C^H?, 

CH^,  H 2 ,  CO,  and  ^0  could  be  adjusted  to  make  the  permeability  increase*” 
or  decrease.  Actual  phenolic  ablation  products,  obtained  from  impregnated 
yams,  were  used  with  the  char  samples.  Their  effect  was  to  decrease  tne 
permeability  of  the  samples.  It  was  also  found  that  the  chars  would  crack 
from  thermal  stresses  when  they  were  exposed  to  high  temperatures  (above 
2800°K)  for  more  than  five  minutes . 

The  analysis  presented  in  the  previous  report  for  the  prediction  of  the 
change  in  permeability  by  carbon  deposition  from  methane  has  been  verified. 
(Reference  16).  Two  figures  are  shown  comparing  this  analysis  with  experi¬ 
ment.  This  analysis  was  also  extended  to  demonstrate  how  it  could  be  used 
to  predict  the  permeability  in  the  char  layer  during  re-entry. 

It  has  been  found  that  the  experiments  with  the  ablative  char  samples  aiouo 
were  not  a  sufficient  basis  on  which  to  base  a  jiecise  analytical  model  of 
carbon  deposition  and  removal  in  ablative  chars.  The  primary  difficulty 
was  the  low  permeability  of  the  chars.  The  low  permeability  required  tne 
use  of  small  gas  flow  rates  which  in  turn  limited  the  experimental  tc.pe na¬ 
tures  to  prevent  premature  reaction  of  the  test  gas.  Even  using  low  tempera¬ 
tures,  premature  reaction  of  the  test  gas  was  frequently  observed  vita  tne 
char  samples.  Using  uniform  more  permeable  porous  carbon  or  graphite 
samples,  the  functional  kinetic  relationships  among  the  various  gases  inter¬ 
acting  with  the  chars  can  be  determined.  Some  results  have  already  been  o  i 
doing  this  and  are  included  in  this  report. 

It  is  recommended  chat  the  analytical  work  which  has  been  initiated  oe 
extended  to  include  increasingly  complex  ablatiun  models,  so  that  maxi¬ 
mum  use  can  be  made  of  the  results  obtained  during  the  present  project. 

In  conjunction  with  this,  further  experimental  data  should  be  obtained  on 
the  kinetics  of  the  reactions  undergone  by  gas  mixtures  typical  of  the  de¬ 
composition  products  of  ablative  materials.  Also  there  is  need  for  more 
detailed  information  on  the  internal  str>actura  of  heat  shield  materials 
at  various  stages  of  decomposition. 


-25/26- 


IX 


REFERENCES 


1.  Weger,  E.,  Brew,  J.,  and  Schvind,  R.  A.,  "An  Investigation  of  Carbon 
Deposition  in  Chars",  BSD  TR  66-385,  Final  Report  USAF  Contract 

AF  05  (694)-695,  Washington  University,  November  1966. 

2.  Darcy,  H.,  "Les  Fontaines  Publiques  de  Ville  de  Dijon",  Dalmont, 

Paris,  1856. 

3.  Scheidegger,  A.  E.,  "The  Physics  of  Flow  through  Porous  Media", 

The  MacMillan  Company,  New  York,  1960. 

4.  Kozeny,  J.,  Akademie  der  Wissenschaften,  Vienna  Almanach,  Mathematisch- 
naturwissenschaftliche  Klosse-Anzeiger-Sitzungsberichte ,  Vol.  Ila, 

No.  136,  (271  -  282),  1927. 

5.  Palmer,  H.  B.,  and  Cullis,  C.  F.,  "The  Formation  of  Carbon  from 
Cases",  Chemistry  and  Physics  of  Carbon,  edited  by  P.  L.  Walker, 

Marcel  Dekker,  New  York,  No.  1,  Chapter  5,  1965. 

6.  Kevorkian,  V.,  Heath,  C.  E.,  and  Boudart,  M.,  "The  Decomposition  of 
Methane  in  Shock  Waves",  The  Journal  of  Physical  Chemistry,  Vol.  64, 
(964  -  968),  August  1960. 

7.  Skinner,  G.  B.,  and  Ruehrvein,  R.  A . ,  "Shock  Tube  Studies  on  the 
Pyrolysis  and  Oxidation  of  Methane",  Journal  of  Physical  Chemistry, 

Vol.  63,  (1736  -  1742),  1959. 

8.  Porter,  G.,  "Carbon  Formation  in  the  Combustion  Wave",  Fourth 
Symposium  or.  Combustion,  (Massachusetts  Institute  of  Technology), 

(228  -  232),  1952. 

9-  Kozlov,  G.  I.,  and  Knorre,  V.  G. ,  "Single-pulse  Shock  Tube  Studies  on 
the  Kinetics  of  the  Thermal  Decomposition  of  Methane",  Combustion 
and  Flame,  Vol.  6,  (253  -  263),  1962. 

10.  Silcocks,  C.  G.,  "The  Kinetics  of  the  Thermal  Polymerization  of 
Acetylene",  Proceedings  of  the  Royal  Society  of  London,  Series  A, 

Vol.  242,  (411  -  429),  1957. 

11.  Kassel,  L.  S.,  "Thermodynamic  Functions  of  Nitrous  Oxide  and  Carbon 
Dioxide",  Journal  of  American  Chemical  Society,  Vol.  56,  (1838  -  1842), 
1934. 


-27- 


12.  Gaydon,  A.  G.,  and  Wolfhara,  H.  G.,  "Flames",  Chapman  and  Hall, 

London,  Chapter  VIII,  1953. 

13.  Sanders,  W.  A.,  "Rates  of  Hydrogen-Graphite  Reaction  Between  1550 
and  2260*C",  National  Aeronautics  and  Space  Administration  NASA  IN 
D-2738. 

14.  Ergun,  S.,  and  Mentser,  M.,  "Reactions  of  Carbon  with  CO2  and  Steam", 
Chemistry  and  Physics  of  Carbon,  edited  by  P.  L.  Walker,  Marcel  Dekker, 
New  York,  Vol.  1,  Chapter  4,  1965. 

15.  Brown,  A.  R.  G.,  and  Watt,  W.,  "The  Preparation  and  Properties  of 
High-Temperature  Pyrolytic  Carbon",  Industrial  Carbon  and  Graphite, 
Papers  Conference,  London  (86  -  100),  1957. 

16.  Schwind,  R.  A.,  "Carbon  Deposition  in  Porous  Media  by  the  Thermal 
Decomposition  of  Methane",  Doctor  of  Science  Dissertation,  Sever 
Institute  of  Washington  University,  January  1968. 


-28- 


TABLE  1 


Equilibrium  Constant  for  Boudouard  Reaction 


T,  (°K) 

Keq 

300 

I.63  x  10"21 

400 

5.40  x  10“l4 

500 

1.82  x  10"9 

600 

1.90  x  10"6 

700 

2.71  x  10‘4 

800 

1.11  x  10'2 

900 

0.195 

1000 

1.91 

1100 

1.23  x  10 

1200 

5.73  x  10 

1250 

1.12  x  102 

1300 

2.09  x  102 

i4oo 

6.26  x  102 

1500 

I.63  x  103 

1750 

1.06  x  104 

2000 

4.21  x  104 

2500 

2.72  x  105 

3000 

8.80  x  105 

Boudouard  React j  on:  COg  +  C 
Equilibrium  Constant:  Ke(}  * 


2C0 


pco2 
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TABLE  2 


Equilibrium  Constants  for  the  Dehydrogenation  Reaction  of  Some  Hydrocarbons 


Hydrocarbon 

Temperature 

1527°  C  (1800°  K) 

2327°  C  (2600°  K) 

Equilibrium 

constant, 

K 

Equilibrium 

constant, 

K 

Methane  (CH.  ) 

4 

7.54X10”4 

1.19X10'4 

Acetylene  (C.^) 

1.98X10”4 

1.81X10"2 

Ethylene  (^H^) 

4.78X10"6 

9.91X10'6 

Methyl  radical  (CH^) 

9.47X10'6 

1.30X10"4 

Methylene  radical  (C^) 

3.50X10"7 

-4 

1.06X10  H 

Methenyl  radical  (CH) 

_ _ _ 

3.67X10"12 

7.08X10”7 
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Data  on  AVCO/RAD  Char rod  Cylinders 
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CO 

supplied  5  July  1967  (DC8  was  an  extra  charred  cylinder) 

approximation  based  upon  Che  downstream  thermocouple  reading  for  DC7,  as  well  as*  the  peak  back  face  temperature 
for  DCS  since  DC7  and  DC8  were  run  at  same  test  conditions. 


AS  RECZIVCD  aiAkACTT*l7.ATl<*i  (CUAAS) 


TABLE  5 

As  Received  Characterization  (Carbons  and  Craphltcs) 
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j. 
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DENSITY  OF  DEPOSIT  GM/CM' 


I 


f 


& 


b. 


TEMPERATURE,  8C 

FIGURE  1:  VARIATION  OF  DENSITY  OF  CARBONACEOUS 
DEPOSIT  WITH  TEMPERATURE 

O  DEPOSITED  FROM  METHANE  AT  15  CM  Hg 
PARTIAL  PRESSURE 

□  DEPOSITED  FROM  BENZENE  AT  2.5  CM  Hg 
PARTIAL  PRESSURE 

A  DEPOSITED  FROM  PROPANE  AT  15  CM  Hg 
PARTIAL  PRESSURE 
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FIGURE  2-  INDUCTION  FURNACE 
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FIGURE  3:  WATER  INJECTION  SYSTEM 


TEST 

GAS 


FIGURE  4:  YARN  APPARATUS 
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OUTSIDE 

SURFACE 


FIGURE  5:  PHOTOMICROGRAPH  OF  R6300-HP-DC5-11 .OFH-XS 

AS  RECEIVED 


DISTANCE  FROM  HOT  END  OF  CYLINDER  (IN) 


FIGURE  6:  COMPARISON  OF  AFTER  CONDITIONING  PERMEABILITIES 


FIGURE  7:  PORE  SIZE  DISTRIBUTION  BEFORE  AND  AFTER  CONDITIONING 


FIGURE  8:  R6300-HP-DC4-0.6FH-XS 

AFTER  CONDITIONING  AT  2875°K 
EXTERIOR  VIEW 
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FIGURE  9:  R6300-HP-DC4-0.6FH-XS 

AFTER  CONDITIONING  AT  2875°K 
INTERIOR  VIEW 


FIGURE  11:  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  MOBILITIES 
IN  PC-60A-I  SAMPLE  FOR  CARBON  DEPOSITION  FROM  METHANE 


1.0 


FIGURE  12:  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  MOBILITIES 
IN  PG-25D-3  SAMPLE  FOR  CARBON  DEPOSITION  FROM  MEfHANE 
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FIGURE  13:  EFFECT  OF  METHANE  ON  R6300-HP-DC5-0 .7FH-XS 


FIGURE  15:  EFFECT  OF  ACETYLENE  ON  R6300-HP-DC1-8.2FH-XS  (MACHINED) 
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FIGURE  16:  EFFECT  OF  ACETYLENE  ON  PG-45D-: 


FIGURE  18:  EFFECT  OF  ACETYLENE  ON  PG-45F-: 


FIGURE  19:  EFFECT  OF  ACETYLENE  ON  PC-45C-2 


OUTSIDE 
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FIGURE  20:  PHOTOMICROGRAPH  OF  R63Q0-HP-DC1-8.2FH-XS 

(MACHINED) 

AFTER  EXPOSURE  TO  ACETYLENE 


FIGURE  21:  PHOTOMICROGRAPH  OF  R6300-HP-DC1-8.2FH-XS 

(MACHINED) 

AFTER  EXPOSURE  TO  ACETYLENE 


FIGURE  22:  EFFECT  OF  WATER  ON’  R6300-HP-DC3-1 .75FH-XS 
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FIGURE  23:  PHOTOMICROGRAPH  OF  R6300-HP-DC3-1.75FH-XS 
AFTER  EXPOSURE  TO  WATER 


AFTER  EXPOSURE  TO  WATER 


FIGURE  25:  EFFECT  OF  METHANE  AND  HYDROGEN  ON  PG-45G-I 


FIGURE  26:  EFFECT  OF  METHANE  AND  HYDROGEN  ON  R6300-HP-DC5-3.55-FH-XS 


METHANE  DATA 


FIGURE  28:  EFFECT  OF  ACETYLENE  AND  HYDROGEN  ON  R6300-HP-DC5-4 .6FH-XS 


FIGURE  29:  EFFECT  OF  ACETYLENE  AND  HYDROGEN  ON  PC-45-G 


FIGURE  31:  EFFECT  OF  ACETYLENE  AND  CARBON 

MONOXIDE  ON  R6300-HP-DC5-8.4FH-XS 
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FIGURE  35:  EFFECT  OF  METHANE,  ACETYLENE,  HYDROGEN,  CARBON  MONOXIDE 
AND  WATER  ON  R6300-HP-DC8-7 .7FH-XS 
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FIGURE  36:  EFFECT  OF  METHANE ,  ACETYLENE,  HYDROGEN,  CARBON  MONOXIDE 
AND  WATER  ON  R6300-HP-DC8-9.9FH-XS 
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APPENDIX:  CHAR  GROWTH  SIMULATION  PROGRAM 


10. I  CARBON  DEPOSITION  EQUATIONS 

The  char  simulation  program  is  based  on  the  solution  of  the  equations 
predicting  the  change  in  permeability  due  to  carbon  deposition  from  methane. 
These  equations  are  listed  below  and  were  derived  in  the  previous  report 
(BSD  TR  66-385) . 


C„  m|-  exp  [C^ 


^  M1/2(rx+y) 


3M  2  *1  Z  ^ 
3t 


dy] 


?  M1/2(r1+y) 

1  +  fx  [1-exp  (CAJ  - q -  dy)] 


(10.1) 


RT 


Q-p  ^  +  n(CH4)l“  n(CH4)1 


(10.2) 


p  2  _  _E1_ 

2  xLK 
o. 


y[nT1  +  U(CH4)1"  n(CH4)3 


(r.  +  y)  M 

l2 


>1/2 


dy, 


(10.3) 


*(CHA)  -  n(CHA)1 


M1/2(rx  +  y) 


4  |  Q 

o 


dy 


(10.4) 


where 


2k  (M.».)c 

n  - - - - 

2  RTp- 


(10.5) 
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The  program  solves  these  equations  for  a  linear  char  growth  rate,  a 
linear  temperature  profile,  and  a  constant  methane  and  inert  gas  flow  rate 
at  a  plane  where  the  gas  enters  the  char  for  a  cylindrical  char.  A  sim¬ 
plified  flow  chart  showing  the  calculation  procedure  is  found  on  page  82 . 

The  input  variables  are  listed  on  page  83 .  The  program  was  written  in 
Fortran  IV  for  use  on  an  IBK  360  computer  operating  with  a  RAGS  system. 

10.2  COMPUTER  PROGRAM  NOMENCLATURE 

The  following  subroutines  and  function  programs  were  used  in  the  char 
simulation  program. 

Function  SUMN(I,  IM) — calculates  the  sum  of  the  squares  of  the  difference 
between  the  last  two  iterations  on  the  number  of  moles  of  methane  at  a 
given  time,  I.  The  summation  is  over  each  radial  position  existing  at  that 
time. 

Function  SUHM (I , IM) — calculates  the  sum  of  the  squares  of  the  difference 
between  the  last  two  iterations  in  the  value  of  the  mobility  at  a  given 
time,  I.  The  summation  is  over  each  radial  position  existing  at  that  time. 

Function  FF(I,J) — calculates  the  right  hand  side  of  equation  10.1  for 
time,  I  and  position  J. 

Subroutine  TEMP — calculates  the  values  of  the  temperature  dependent 
variables  for  each  radial  position  existing  at  time  I. 

Subroutine  PC2 — performs  a  second  order  predictor-corrector  solution  of 
equation  10.1.  Also  includes  Euler's  method  followed  by  second  order 
corrector. 

Subroutine  PC4 — performs  a  Hamming's  fourth  order  predictor-corrector 
solution  of  equation  10.1. 

Subroutine  RFND — calculates  the  root  of  equation  10.4  for  given  values 
of  mobility  and  pressure.  The  secant  method  is  used  to  find  the  root. 

Subroutine  1NTEG — performs  integration  in  equation  10.4  by  trapezoidal 
rule. 

Subroutine  PRES— solves  equation  10.3  for  given  values  of  the  moles  of 
methane  and  the  mobility  for  each  radial  position  at  time  I.  Again  the 
trapezoidal  rule  was  used. 
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10.3  NONSUBSCRIPTED  VARIABLES 


ANAR  molar  flow  rate  of  inert  gas 

ANCH4  molar  flow  rate  of  methane  gas 

AVM  average  mobility  at  any  time 

CCI  convergence  criterion 

CC2  square  of  convergence  criterion 

CSTEP  determines  whether  step  size  is  to  be  changed 

C2  constant  defined  by  equation  10.5 

C4  constant  defined  by  equation  10.6 

DELT  time  step  size 

DELY  step  size  in  y  direction 

DTDXH  constant  temperature  gradient  across  the  char 
DTMX  final  step  sire,  maximum  step  size 

HEAD  sample  designation 

HI  height  of  the  sample 

I  index  to  keep  track  of  times  stored  in  the  computer 

IM  index  to  keep  track  of  actual  time  that  the  times  stored  in 

the  computer  represent 

JS  counter  to  keep  track  of  the  value  of  the  index  corresponding 

to  the  inside  surface  of  the  char 

MPOUT  determines  whether  intermediate  pressure  and  number  of  moles 
of  CH^  output  is  printed  out 

KSEP  determines  if  mobilities  and  pressures  ac  each  y  should  be 

printed  out  at  selected  increments  In  time 

allows  mobilities  at  each  y  for  all  times  to  be  printed  out 


NR 


input  tape  number 
NS  indicates  new  sample  follows 

NTH  indicates  the  number  of  increments  in  y  which  are  used 

NW  output  tape  number 

N2  indicates  type  of  inert  gas  used 

ODIA  outside  diameter  of  the  sample 

PERO  initial  permeability  of  the  sample 

PORO  initial  porosity  of  the  sample 

P2B  atmospheric  pressure 

KABIF  inside  radius  of  the  sample  in  feet 

RK  rate  constant  for  the  decomposition  of  methane 

ROAR  inert  gas  flow  rate 

R0CH4  methane  gas  flow  rate 

T  temperature  of  the  sample  during  the  carbon  deposition 

TH  thickness  of  the  sample 

THEND  maximum  thickness  of  the  char 

TIME  total  time 

TMAX  maximum  time  of  the  deposition  experiment 

THPO  time  increment  for  normal  output 

TREF  temperature  of  outside  surface  of  the  char 

TREFX  temperature  of  inside  surface  of  the  char 

VISAR  viscosity  of  the  inert  gas 

VISCH4  viscosity  of  the  methane 

viacoslty  of  the  inert  gas-methane  mixture 


VISMX 


10.4  SUBSCRIPTED  VARIABLES 

BM(I,  J)  mobility  at  time  I  and  position  J 

BNCH4 (I ,  J)  number  of  moles  of  methane  at  time  I  and  position  J 

FNCH4(J)  number  of  moles  of  methane  at  position  J 

P(I,  J)  pressure  at  time  I  and  position  J 

RK(I,  K)  rate  constant  for  decomposition  of  CH,  at  time  I  and 

position  K 

SINT  (I,  J)  value  of  integral  in  equation  10.4  for  time  I  and 

position  J 

SPINT(K)  integral  of  equation  10.3 

T(I,  K)  temperature  of  char  at  time  I  and  position  K 

VISKX(I,  K)  viscosity  of  the  methane-inert  gas  mixture  at  time  I 
aud  position  K 


I 


FIGURE  A-l:  SIMPLIFIED  FLOW  CL. ART  FOR  CHAR  SIMULATION  PROGRAM 
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INPUT  FOR  CHAR  SIMULATION  COMPUTER  PROGRAM 


41-42  X  NTH  Number  of  Inerii.ijic#  used  to  reach  maxlnum  char 

thickness  (oust  be  an  even  number'/ 


COMPUTER  LISTING  OF  CHAR  SIMULATION  PROGRAM 


C  CHAR  SIMULATION  PROGRAM 

C  EXACT  THEORY  T  =  T (TIME/ Y)  THICKNESS  VARIES 

C  THIS  PROGRAM  SOLVES  FOR  TliE  M03ILITY,  PRESSURE, AND 

C  THE  MOLES  OF  CH4  IN  THE  CHAR  OF  A  REENTERING  ABLATOR 

C  WHERE  THE  CHAR  THICKNESS  IS  ALLOWED  TO  VARY  FROM  0  TU  TH. 

C  INPUT  IS  READ  INTO  THE  COMPUTER  AS  FOLLOWS 
C  CARO  1.  COL.  1-60  HEADING,  A  FIELD 

C  CARD  2.  COL.  1-2  INERT  GAS  USED  (0  FUR  AkGON,  1  FOR  HE, 

C  2  FuR  H2 

C  CARD  3.  COL.  1-10  Cli4  FLOW  RATE,  COL.  11-20  INERT  GAS 

C  FLOW  RATE,  COL.  21-30  OUTSIDE  DIAMETCR,  CUL.  31-40 

C  HEIGHT,  COL.  41-50  FINAL  THICKNESS  OF  THE  CHAR, 

C  FORMAT ( 5D10 . 4 ) 

C  CARD  4.  COL.  1-10  LENGTH  OF  TIME  FOR  THE  THICKNESS 
C  TO  REACH  A  STEADY  STATE  VALUE 

C  CARD  5.  COL.  1-10  MAXIMUM  TEMPERATURE  THAT  TliE 
C  OUTSIDE  SURFACE  OF  THE  CHAR  WILL  REACH,  CUL.  11-20 

C  TEMPERATURE  OF  INSIDE  SURFACE  OF  THE  CHAR  (USUALLY 

C  TAKEN  AS  THE  DECOMPOSITION  TEMPERATURE  OF  PHENOLIC  RESIN), 

C  COL.  21-3C  PERMEABILITY  OF  CHAR  WHICH  IS  ORIGII  ALLY 
C  FORMED,  COL.  31-40  OR  I  GIAN L  POROSITY  OF  THE  CHAR, 

C  COL.  41-42  HUMBER  OF  STEPS  USED  ACROSS  CHAR 

DOUBLE  PRECISION  P(  5,  21 ) ,  SNCH4  (  5  ,21 )  ,BM(  5  ,21 )  ,AFF  ( 21)  , 

1  SIM  T(  5,21)  ,TI5,21),RK(5,21  ),VI  SMX(5,21)  ,FNCH4(21), 

2GM{  21 )  ,ROCH4^ 

3001  At  HT,  THENO,  OELT  ,XCH4»RA01»RAD1F, VCH4M,  ANCH4,  VN2M  , 

4ANA  R,  VARM, 

5VHE  Mt  ANTOT  ,  AN'TCH ,  PERQ,  PORO , TREF  ,DTDTH *DELY»TMAX»  ANTH, 

6C5  ,C6  ,ROAR 
7 , XAR , T  IME  ,PREO ,TREF I 
DIMENSION  HEAD(  15  ) 

COMMON  P  ,BNCH4 , DM ,A FT , SINT , T » RK » VI SMX »  FNCH4  ,  GM ,  C6  ,C5 , 
1XCH4,NTH 
NR=  l 
NW=  2 

500  REAO(  NR, 5  )  HEAD 
5  FORMAT!  15A4) 

C  N2=  0  FOR  ARGON  =1  FOR  HE.  *2  FOR  N2 
READ  (MR,  105)  N2 
1C5  FORMAT (12) 

N  2=  N  2  + 1 

WRITE  ( NW , 40)  HEAO 
40  FORMAT ( 1H1, 15A4, ///) 

C  SAMPLE  DIMENSIONS  ARE  IN  l.'CPES 
C  FLOW  RATES  ARE  FT3/MIH 

READ  (1,10)  RQCH4yR0AR,0DIA,HT, THENO 
10  FORMAT (5010.4) 

C  TMAXslO  ,20,  30  SEC. 
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C  TMA X  =  LENGTH  OF  TIME  FOR  THICKNESS  TO  REACH  STEADY 
1STATE  VALUE 
READ  (1,16)  TMAX 
16  FOR MAT(D 10.4) 

RAD  1=  .  500*00 IA-THEND 
RAD1F=RAD1/12. 

VC!!4n  =  R0CII4 

AKCH4  =  VC!.‘4M  *  .  9025G2G000/60 . 

GO  TO  (C/12/9)/N2 
9  VN2I1  =  ROAR 

ANAR  =  VN2M  *  (7,20350-2)/ (28.  *60.) 

GO  TO  11 
8  VARH  =  ROAR 

ANAR  =  VARH  *  ( 2 . 554D-3 )/6 0 . 

GO  TO  11 
12  VI ! Eii  =  ROAR 

ANAR  =VHEM*(  2.55D-3)/60. 

U  ANT0T=A.NCH4+ANAR 

WRITE  ( NW  ,45  )  ANCH4  ,ANAR 

45  FORMAT  (8H0NCH4  =  ,D14.6,8HPER  SEC  ,f'tX,6HNAR  -  ,Dl^6, 

18H>ER  SEC.,/) 

ANTCH=ANT0T+ANCH4 

C  NTH=  NUMBER  OF  INTERVALS  IN  Y-DIRECTION  MUST  36  EVEN 
1N0. 

C  TRflF  IS  TEMPERATURE  OF  OUTSIDE  SURFACE 
C  Titr.FI  IS  TEMPERATURE  OF  PNENOLIC  DECOMPOSITION 
1  READ  (NR,  15)  TREF,TR£FI ,P£RO,PQRO,NTK 
15  FORMAT (4D 10. 4,12) 

XCH4=ANCH4/ANT0T 
XAR  =  ANAR/ANTOT 
AMTH=DFLOAT ( NTH ) 

OTOTH=  (  TR EF-TREF I  )/THEN0*l2. 

DELT=TMAX/ANTH 
OELY=THENO/( 12.*ANTH) 

WRITE  ( NW ,30 )  TREF , THEND,  TMAX 
30  FORMAT ( 9H-TEMP .  =  ,DL0.4.5X,7H  TH.  »  • DIO. 4 .AH  TMAX  » 
If  01  0.4,//) 

C5=-24.*2000./140.*XCH4/39.35D0 
C6=-6.28  3185307*HT*P0R0/12. 

NTH=NTH+1 
I  M=  1 
1  =  1 

BM(  1,NTH)  =  1. 

8NC  H4 ( 1 ,NTH  )=ANCH4 
S  INT  (  1,NTH)  =  0. 

FNC  H4 (NTH )* ANCH4 

J  =N  TH 

TIME-0. 

CALL  TEMPI  I,  IM,  TIME, DELYfDTDTHt  ANCH4»ANT0T,  ANAR, XARt 
1HT,N2) 

Call  PRES  (DELY,  ANTCHfRADIFfI  flM.HTfPERO) 


-85- 


WRITE  (NW,68)  T  IM6  ,BM{  I  M,  NTH )  ,  P  ( I  H,  NTH)  ,  BNCH4  ( IM,NTH ) 
68  FORMAT (8H  TIME  *  , 010 .4/ 1  3020 .8  )> 

I  M=2 
1=2 

RAD  IF=RAD  1F-DELY 
JS=NTH-1  M+l 
BM(  I,JSi=l. 

BNCH4I I f  J  S )=ANC  H4 
FNCH4  { J  S)=ANCH.4 
TIME=TIM£+DElT/60. 

CALL  TEMPII , IM, TIME, DELY, DTDTH,ANCH4 ,ANTQT, ANAR, XAR„ 
1HT,N2) 

J  =M  FH 

pre  n=i. 

CALL  PC2  II,IM,J,?RED,ANT0T,ANCH4,RAD1F,DELY,  ANTCH,HT, 
1PER  0»DE  LT  ) 

GO-  TO  151 

111  IF  (SUMMI I»IM)-C i.OD-26))  151,151,152 

151  J  S=NTH-  IM  +  1 

WRITE  (NW,68)  T  IME,  (  BM(  I » JJ) »  PI  I ,  JJ  ) ,  BNCH4 1 1,  JJ  )  , 

1J  J=  J  S,  NTH ) 

1=1+1 
IM= I M+  1 

RAO  1F=RAD1F-DELY 
IF  l  1-6)  108,105,109 
108  J  S=NTH- 1 M+ 1 
BM( I , JS)= l • 

GM{ J  S)  =  l • 

BNCH4 ( I , JS)=ANCH4 
FNC  H4 1 JS  )=ANCH4 
KK= JS+l 

DO  201  K=KK ,  NTH 
BMl  I , K )=BM( I— 1, K ) 

2C1  BNCH4 I I , K )=BNCH4( I- 1 , K) 

IF  { I M-NTH)  116,116,112 
lit  T  flJC-T  T  0ELT/60  * 

CALL  TEMPII ,IM, TIME, DELY, DTDTH,ANCH4,ANT0T, ANAR, XAR, 
1HT,N2) 

152  J  S=NTH-  IM+2 

DO  211  K=JS,NTH 
211  GM( K )  =  8M I  I , K 1 
J=J  S 
PR6D=1. 

CALL  PC 21  I  ,U-!,J»PRED,ANTOT,ANCH4,RAOIF,DELY,ANTCH,HT, 
IP  ERO  ,DELT ) 

J=J  +  1 

1  C6  PRED=2* 

CALL  PC2  1 1 , 1 M, J , PRED, ANTOT « ANCH4,RAD1F,DELY,  ANTCH,HT  , 
1PER0, DELT ) 

IF  IJ-NTH)  101,111,111 
101  J=J  +  i 


-86- 


IF  { J- { NTH- IM+4 I )  106,106,107 
1C7  CALL  PC4I I ,IM»J,PRED»ANTQT,ANCH4,RAC1F»D£LY*ANTCH»HT» 
1PJS0  , DELT) 

IF  (J-NTH)  101,111,111 
109  1=0 

JS=  NTH-IM+6 
00  200  K=2» 5 
1=1  +  1 
JS=JS-1 

00  200  J=JS,N7H 
BMI  I  ,.J )  =  BM»  K ,  J ) 

P(I  » J ) =P I K, J  ) 

SINK  1 , J)=SINT{K* J) 

T ( I , J)=T IK,  J) 

RKI  I » J )=RK{  K,J) 

VISMXI  I ,  J  )=VISMX(  K,  J) 

2C0  BNCH4(I,J)=BNCH4(K»J) 

1=5 

GO  TO  108 
112  GO  TO  500 
600  END 
/FTC 

DOUBLE  PRECISION  FUNCTION  SUMN(I,IM) 

DOUBLE  PRECISION  P(  5  ,2  1 1  »BNCH4t  5, 21),  BM  (5, 21),  AFT  (21), 
IS  INK  5,21),TI5,21),RK(5,2i),VI$MXI5,21),FNCH4(21), 

2GM(  21) 

3  ,C6  ,C5,XCH4,X,Q 

COMMON  P  ,BNCH4»BM»AFT , SINT,T,RK,  VI  SMX,  FNCHA  ,GM,C6  ,C5, 
1XCH4»NTH 
X=0. 

JSS  =  NTH-I  M+2 
DO  61  K=JSS , NTH 
IF  (BNCH4I I  ,K) )  63,63,62 
63  Q=0 . 

GO  TO  61 

62  0=(  FNCH4  (K  )  — 8NCH4I  1  ,K)  )/BNCH4(I»K) 

61  X=X+Q*Q 
SUMN=X 
RETURN 
END 

/FTC 

DOUBLE  PRECISION  FUNCTION  SUMMII,IM) 

OOUBLE  PRECISION  PI  5 ,2 1 ) ,BNCH4(  5, 21 ) ,  BM (5, 21) ,  AFT! 2 1 ) » 
IS  INTI  5,21)  ,T(  5,21)  ,RK{  5, 21),  VISMXI  5,21)  ,FNCH4  (21)  , 

2GMI 21), 

3C6,  C5 » XCH4,  R,  Y 

COMMON  P,BNCH4, BM, AFT, SINT, T,RK , VISMX ,FNCH4 ,GM, C6 ,C 5, 
1XCH  4 ,NTH 
Y=0  . 

J  SS=NTH-  IM+2 
DO  81  K*J  SS,NTH 


« 
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R=(  GMI  K)— SMI  I ,K) J/BMII ,K) 

81  Y=Y+R*R 
SUMM=Y 
RETURN 
END 

/FTC  REF 

DOUBLE  PRECISION  FUNCTION  FF(I,J) 

DOUBLE  PRECISION  PI  5,  21) » BNCH4  { 5,  21 )  ,  BH  15 ,2  1) ,  AFTI  21 )  » 
1SINT(5,21) ,TI5,21),RK(5,21),VISKX(5,2l) ,FNCH4(21), 

2GM( 21 ) 

3,C6,C5,XCH4,  Z 

COMMON  P,BNCH4,  BM,  AFT,  SINT, T»RK  ,VISMX  ,FNCH4,GM,C6,CS 
1XCH4.NTH 

Z=D  EXP  I  C6*S  INTI  I,  J)  ) 

FF=C5#8M(I  ,J)*P(  I,J)*RKI  I,J)*Z/{T(:»J)*TII,J)* 

II  1.  +  XCH4*(1.-Z)  )) 

RETURN 

END 

/FTC 

SUBROUTINE  TEMPI  I ,  I  N,  TI  ME  ,DELY  t  OTDTH,  ANCH4 ,  ANTOT,  AN  AR , 
1 XAR  , HT,N2  ) 

DOUBLE  PRECISION  PI  5 ,21 )  ,BNCH4(  5,21 )  ,  BM  (5 ,2 1 )  ,  AFT  (21 )  ♦ 
IS  INTI5,21),T(5,2l),RKI5,21),VISMX{5,21) ,FNCH4(21), 

2GMI  21)  ,C6,C5, 

3DTD  TH.ANCH4 , ANT OT  f  ANAR, HT,Y,  VI  SCH4 ,  VISAR  ,XAR  ,PH  12, 

4PH2  I ,  XCH4  » TIME 
5,  DEL Y 

COMMON  P,  BNCH4,  BM  »AFT ,SINT ,T  ,RK,VISMX,FNCH4,GM,  C6,C5, 
IX  CH  4 ,  NT  H 
K=NTH 
IY= IM- 1 

12  Y=DELY*DFLOATIIY) 

TII  ,K)=DTDTH*Y+2.03 
VISCH4=  1 2. 930-2)  *(6.71970-4) 

C  THIS  IS  VISCOSITY  OF  CH4  AT  1000  K  (PERRYS)  LBM/ 

IF T- SEC 

GO  TO  I  212,216,213) ,N2 

213  V  IS  AR  =  I  6.7 19  7D-8)  *174.  *397.  *  i  (T  ( I,  K  )/293.  )**1 .5 )/  (T  (I , 
IKK  104.) 

GO  TO  214 

212  VISAR=(6.7170-8)*22  2.*435  .*l(T(  I  ,  K )/ 293. )  **1.  5)  /  ITU, 
1KJ+142. ) 

GO  TO  214 

216  VISAR -=(6.7 1970-8  i  *1 94.  *363.  *  (  (T  (I  »K)  /293.  )  **l»5  )/  (T  I  I , 
1KK70.) 

C  1=CH4  AND  2=  N2  OR  AR 

214  PHI 1= 1. 

PH l 2« .353550 C/11.4018 00 **.5)*( !.♦( 2. 48900**. 25) * 

II  IV  ISCH4/VI  SAR) 

2**.S) )**2. 
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PH2 1=. 3535500/(3. 489D0**.5)*{ l.+C .401800**. 25)* 

1( (V ISAR/V1SCH4) 

2**. 5) )**2. 

PH  22=1. 

VI SMX ( I , K }»XCH4*VISCH4/ (XCH4* PHI 1+XAR*PH12 ) *XAR  *V IS AR/ 
11  )CH4*PH?l 
2+XA R*PH22) 

RK(  I,K)*10.**(15.160C«2.28D4/T(  I,K> » 

K=<—  1 
I  Y=  IY—  1 

IF  (IY)  11,12,12 
11  RETURN 
END 

/FTC  REF 

SUBROUTINE  PC2(  I ,  I  K,  J ,  PRED,  ANTO  T ,  ANCH4  ,RADl  F,DELY  , 

1 ANTCH»HT ,PERO , 

20aT) 

DOUBLE  PRECISION  P(  5 ,2 1 )  ,8NCH4(  5,  21 ) ,  BK  (5,2 1 )  ,  AFTI21)  , 
1SINT(5,21)  ,T(5,21  J  ,RK(5,21>  ,VISPX(5»21)  ,FNCH4(21>, 

2GK(  21),C6,C5, 

3PN1  ,PRED,ANTCT,ANCH4,RAOIF,OELY,ANTCH,HT,PER0,XCH4, 

4  D  FL  T  ,  FM 

COMMON  P,BNCH4,BM,AFT,SINT,7,RK,VISMX,  FNCH4  ,GM,C6  ,C5  , 
1XCH4,  NTH 

IF  (PRED-l.l)  51,51 v6l 

C  FOLLOWING  USES  EULERS  METHOD  TO  APPROX.  M  AT  T*T+DELT 
51  8NC  H4(  I  ,  J}=8  NCH  4  (I  — 1 ,  J) 

TRUN=0. 

69  8M(  I,  J  }  =  8M{  1-1,  J)*DELT*FF  (I-l,J) 

GO  TO  71 

C  NOW  USE  2  NO.  ORDER  PREO. 

61  BNCH4(I  ,J)  =  8NCH4(I-i,J) 

TRU  N=  1. 

BM{  I » J  )  =B  M  (  1-2, J)+2.*DELT*FF  (I-l,J> 

PN1  =  BM(  I,  J) 

71  CALL  PRES  (OELY  ,ANTCH,RAD1F,  I » IM,HT,  PER.Q) 

CALL  RFNOIDELY,  I,  IM  ,ANTOT  ,ANCH4,RA01F) 

IF  (SUMNd,  IM  )  -  (1 .  OC—26)  >  74,71,71 
C  NOW  USE  2  NO  ORDER  CORRECTOR 

74  FM=BM( I- 1 , J ) +DE  LT/2.*(FF(I-l,J)+FF( I , J) ) 

IF  { OABS(  FM-BM(  I , J) )  — 1 1. OD— 13) )  80,73,73 
73  8M  (  I ,  J  )  =  FM 
GC  TO  74 
60  BM{  I ,  J )  =  FM 
85  IF  (TRUN)  681,685,681 
681  BM(  I,  -J)  =B  M(  I ,  J  )  + 1.  /5.  *  (PN1-8M(  I,J)  ) 

685  CALL  PR  ES(  DELY  ,ANTCH » RAD1F,  I,IM,HT,P  ERO) 

CALL  RFNOIDELY,  I,  I M,ANT0T,ANCH4,  RADI  F> 

IF(  SUMNd,  IM)-(  1.00-261)  87,685,685 
87  RETURN 
END 
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FTC  REF 

SUBROUTINE  PC41  I , IM,  J,  FRED , ANTOT,  ANCH4,  RADI  F,  DELY  , 

1 ANT  CH » HT»  P  ERO » 

20ELT) 

DOUBLE  PRECISION  P( 5,21 ) , BNCH4 (5  ,21 ) ,BM< 5,21 ) ,AFT < 21 )  , 
1RK{ 5,21), 

2SINT{5,21),TI5,21),VISKX(  5,21)  ,FNCH4  (21) ,  GM  (21)  ,X  ,C6  , 
3C5,  XCH4 

4, ANT CT,  ANCH4,  RADI F,  DELY  ,ANTCH,HT ,  PE RQ,  PN1,C,PRES»  CELT 
COMMON  P, BNCH4,  8M,  AFT,  SINT,  T,RK  ,VISMX  ,F  NCH4,GH,C6,C5, 
1XCH4  ,NTH 
K004=0 

BNCH4I  I,  J  )=BNCHA(  I—  1 ,  J ) 

BM ( I ,  J )  =  8M  (  1—4  » J  )  *4  *  *OE  LT/  3 •* (  2.*FF  (I -1,  JJ-FFU-2  , 

1J  )+  2.*FF  ( I— 3 ,  J)  ) 

PN1=8M (I ,  J). 

IF  { MOD 4)  52,392,92 
392  X=0. 

92  BK( I,J)=PN1-112./121.*X 

51  CALL  PRES(DELY»ANTCH»RAD1F, I , IM,HT, PERO) 

CALL  RFNDIDELY ,  I  #  I M, ANTOT  ,ANCH4, RADIF) 

IF  { SUMN ( I , I M}— (1 «OD— 26) )  78,51,51 
78  C4TM=0. 

178  C=l./8.*(9.*BMU-l,J)-BMlI-3,J)+3.*DELT*(FFn,J)+2.* 
1FF{  I— 1, J) 

2*-FF  (  I~2,  J  ) )  ) 

IF(  DABSt  C-BMl  I,J)  )-{  1.00-13))  170,76,76 

76  3M(I,J)=C 
C4TM=C4TM+1. 

IF  (C4TM-10)  178,178,602 
602  WRITE  (NW,610)  C,GM(J) 

610  FORM  AT (47H  4  TH*  CRDER  CORR.  DOESN’T  CONVERGE  IN  10 
1CYCLES, 2015*3} 

170  BMl  I ,  J)  =  C 

82  CALL  PRES  (DELY,  ANTCH,  RADI  F,  I ,  IM  ,HT ,PERO ) 

CALL  RFND(OELY,I , IM, ANTOT, ANCH4 , RADIF ) 

IF  (SUMNt  I»IM)-t  l.OD-26))  77,82,82 

77  X=PN1— C 

BM(I,J)  =  C+9./121.*X 

81  CALL  PRES  (DELY,  ANTCH, RADI  F,  I  »  IM,HT,  PERO) 

CALL  RFNDIDELY, I  ,IM,  ANTOT, ANCH4  , RADIF) 

IF  (SUMNI I»IM)— (1*00-26))  195,81,81 
195  REfURN 
END 

/FTC 

SUBROUTINE  RFND  ( DELY ,  I  ,  I M,  ANTOT  ,  ANCH4,  RADIF) 

DOUBLE  PRECISION  PC  5,211 ,  BNCH4J  5,21)  ,BM(  5, 21),  AFT  (21), 
1SINT(5,21),T(5,21),RK(5,21)  ,  VI  SMX(  5  ,21 )  ,FNCH4l  21) , 

2GM(  21 ) ,FNC ( 2 ) , 

3  ANT  CH,  0  EL  Y,  ANTOT  ,ANCH4,  RAD  IF ,  Y,  CNCH4,  DNCH4,  C6 ,  C5  ,XC  H4 
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COMMON  P,BNCH4,BM,AFT,$INT,T,RK,VISMX,FNCH4,GM,C6,C5, 
1XCH4,NTH 
ANTCH-ANT0T+ANCH4 
CON V=0. 

NC0N=0 
JSR=NTH-IM+1 
Y=0  . 

AFT( JSR)=39.35DC*T(I,JSR)/P(I,JSR)*ANT0T 
SINTl  I ,  JSR  )=0  • 

IF  (1-1)  201,201,211 

211  J  $R=J  SR+1 

DO  200  L=JSR ,NTH 
I RT  =0 
Y=Y  +DELY 
CNC  H4=B  NC  H4 ( I , L ) 

203  IRT  =  IRT  +  1 

GO  TO  (205,204) , IRT 

205  AFT (L )=39.3500*T( I ,L )/?( I ,L) *(ANTCH-CNCH4> 

CALL  INTEGII ,L, OELY ,Y,RAD  IF ) 

GO  TO  (208,204) , IRT 
2C8  CNCH4=ANCH4*DEXP(C6*SINT(I ,L) ) 

FNC  (  1  )=8NCH4(  I ,  LJ-CNCH4 
GO  TO  203 

204  FNC (2)=CNCH4-ANCH4*0EXP(C6*SINT ( I ,L) ) 

IF  (FNC(2)-FNC( 1) )  212,213,212 

212  DNC H4=CNCH4— (CNCH4— BNCH4 ( I ,L) ) / ( FNC { 2 )— FNC ( 1 ) )*FNC ( 2 ) 
CON  V=CONV+ 1 . 

IF  ( CNCH4  )  305,207,305 

305  I F  (  OABS  (  (  CNCH4— DNCH4 )  /CNCH4 )  —  ( 1*00—13)  )  2  07  ,  2  06,206 

206  FNC  ( 1  )  =  FNC ( 2 ) 

B(CH4(I,L)  =  CNCH4 
CNC  H4=D  NCH4 

I  F(  CONV-  30  •  )  20  5,  20  5,209 

209  NCQN=NC0N  +  1 

WRITE  ( N  W,2 10  )  CNCH4  ,BNCH4(  I  ,  L) 

210  FORMAT  (2CH  ROOT  NOT  CONVERGING,  2015.  8) 

IF(  NC ON-4)  2  07,207,200 

213  DNC  H4=CNCH4 
2  07  OCN  V=0  . 

BNC  H4  (  I  ,L  '=DNCH  A 
2C0  CONTINUE 
201  RETURN 
END 

/FTC 

SUBROUTINE  I  NTEG  ( I  ,  L,  DEL  Y,  Y,  RADI  F) 

COUBLE  PRECISION  P(  5,21 )  ,  BNCH4  (5 ,21 )  ,  BM  (5  ,2  1  )  ,AFT  (  2  l ) , 
1SIN  T( 5  ,21 )  ,T(5, 21)  ,RK(  5,21 ) ,DELY, Y, RADI  F,  XR AD, DINT  ,A  ,B 
COMMON  P,  BNCH4,BM, AFT, SINT  »T,RK 
XRAD=RAD1F+ Y 

A  *3  SQRT  {  BM(  I  » L-l )  )*  (XR  AD-DEL  Y)/AFT(  L- 1 )  *RK(  I  ♦  L-l ) 
ft -0  50*7(8  Ml  I  ,LI  )  *  XRAD/AFT  (L )  *RK  (  I  ,L  ) 
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DINT=0£LY/2.*(A+8J 
SINT(I,L)*SINT( I,L-1)+0INT 
RETURN 
END 

/FTC 

SUBROUTINE  PRES (OELY, ANTCH, RADI F, I , IM ,HTtPERO I 
DOUBLE  PRECISION  Pf  5,2  I > ,  BNCH4C  5, 21  ) ,  RM  (5  ,21 )  ,AFT(  2  i )  t 
IS  INK  5,21)  ,T(5,21)  ,RK{5,21  J  ,VISMX  (  5,21)  »GM(2 1 ) , C6,C5, 
2XCH4,FNCH4{21)t  SP IN T( 21 ) »DELY  ,ANTCH  » RAD IF  *HT»  PERO* 
3XRAD  »AfB  »OPI  NT* 

4ANTHf Y 

COMMON  P  » BNCH4  »BK»AFT »S  INT»T»  RK»  VI  SMX»  FNCK4  »  GM,  C6  »C5» 
1XCH4*  NTH 

FNCH4 (NTH)=BNCH4( I, NTH) 

ANTH=OFLOAKNTH) 

Y=D  ELY*ANTH 
Pt I ,NTH)=29* 52100 
IF  (1-1)  410» 41C»405 
405  SPINT (NTH)=0* 

K=N  TH-1 
4C0  Y=Y-DELY 

XRAD=RAD1F+Y 

A=( AN  TCH-BNCH4I I « K+ 1 ) )  /  { XRAD* BM  U ,K  U ) ) *V IS  MX ( I , K+l ) * 

1 T( I »K  +  1 ) 

B=(ANTCH-BNCH4(I  ,K)  )/( (  XRAD-OELY)*BM  (I » K) )  *VISMX(  I,K)* 
IT  (I  i  K ) 

DPI  NT=DELY/2.*(A+B) 

SPI  NT  {  K)=S  PI  NT  (K+ 1)  +  0P INT 

?( I » K )=DSQRT ( £95. 26600+ (6 *601640- 2)/ ( HT  *PER  0) * 
lSPINT(K)) 

FNCH41 K )=BNCH4( I »  K) 

IF  (K—(NTH“  IM  +  i  ) )  410*410,420 
420  K=K-l 

GO  TO  400 
410  RETURN 
END 
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